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ABSTRACT 

We use a suite of semi-empirical models to predict the galaxy-galaxy merger rate and relative contributions 
to bulge growth as a function of mass (both halo and stellar), redshift, and mass ratio. The models use empirical 
constraints on the halo occupation distribution, evolved forward in time, to robustly identify where and when 
galaxy mergers occur Together with the results of high-resolution merger simulations, this allows us to quantify 
the relative contributions of mergers with different properties (e.g. mass ratios, gas fractions, redshifts) to the 
bulge population. We compare with observational constraints, and find good agreement. We also provide 
useful fitting functions and make public a cod43 to reproduce the predicted merger rates and contributions to 
bulge mass growth. We identify several robust conclusions. (1) Major mergers dominate the formation and 
assembly of ^ L» bulges and the total spheroid mass density, but minor mergers contribute a non-negligible 
^ 30%. (2) This is mass-dependent: bulge formation and assembly is dominated by more minor mergers in 
lower-mass systems. In higher-mass systems, most bulges originally form in major mergers near ^ L^,, but 
assemble in increasingly minor mergers. (3) The minor/major contribution is also morphology-dependent: 
higher B/T systems preferentially form in more major mergers, with B/T roughly tracing the mass ratio of 
the largest recent merger; lower B/T systems preferentially form in situ from minor minors. (4) Low-mass 
galaxies, being gas-rich, require more mergers to reach the same B/T as high-mass systems. Gas-richness 
dramatically suppresses the absolute efficiency of bulge formation, but does not strongly influence the relative 
contribution of major versus minor mergers. (5) Absolute merger rates at fixed mass ratio increase with galaxy 
mass. (6) Predicted merger rates agree well with those observed in pair and morphology-selected samples, 
but there is evidence that some morphology-selected samples include contamination from minor mergers. (7) 
Predicted rates also agree with the integrated growth in bulge mass density with cosmic time, but with factor 
^ 2 uncertainty in both - up to half the bulge mass density could come from non-merger processes. We 
systematically vary the model assumptions, totaling 10^ model permutations, and quantify the resulting 
uncertainties. Our conclusions regarding the importance of different mergers for bulge formation are very 
robust to these changes. The absolute predicted merger rates are systematically uncertain at the factor ~ 2 
level; uncertainties grow at the lowest masses and high redshifts. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: active — cosmology: theory 



1. INTRODUCTION 

In the now estab lished ACDM cosmology, structure grows 
hierarchically (e.g. I White & Rees 1978), making mergers an 
inescapable element in galaxy formation. Thirty years ago, 
l^omre (1977) proposed the "merger hypothesis," that ma- 
jor mergers between spirals could result in elliptical galax- 
ies, and the combination of detai led ob servations of re- 
cent merger remn ants (ISchweize rl Il982t [L ake & Dresslei 
1986'; 'Dovon et al."1994; Shier & Fischerl998; James et al 
1999; Genze letal. 2001; Tacconi et al. 2002; Dasvraet^ 
20061 12007t iRothberg & JosephI 12004 I2006at Ivan Dokkum 



120051) and e.g. faint shells and tida l features around 
ellipticals jM aUn & Carteil fT980l fT983l i Schweizeii fT980l 
ISchweizer & Seitzer 1992; Schweizer 1996) h ave lent consid - 
erable support to this picture (e.g. B arnes & Hernquistlll992h . 

Mergers are also linked to starburst galaxies and lu- 
minous quasars. By exciting tidal torques that lead 
to rapid inflows o f gas into the centers of galaxies 
(iBarnes & Hernquistl[l99U 11996), mer gers provide t he fuel 
to power intense starbursts (Mihos & Her nquisfl I l994t 



1996) 
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to feed rapid black hole growth (Jbi Matte o et al.l 
iHopkins et al.l l2005bllah . and, through various asso- 
ciated feedback channels, convert blue, star-f orming galax- 
ies into quiescent, red ones (e.g. ISpringel et al.l 1200 5a b). 
Mergers are inevitably associated with the most lumi- 
nous star-for ming systems, fro m ULIRGs in the lo- 
cal Universe (ISoifer et al l Il984allbl; [Joseph & Wright [T98l 
ISanders & Mirabel! Il996h to bright sub-mi l limeter galax- 
ies at high redshifts (Alexander et al. 2005'; Yo unger et al.l 
2008b; Shapiro et al. 2008; Tacconi et al. 2008), and proper- 
ties ranging from their observed kinematics, structural cor- 
relations, and cl ustering link thes e populations to massive 
elUpticals today jLake & Dressle3 [l986; Dovon etal. 1994; 
IShier & Fischer '1998'; 'James et al.^ 1999; Genzel etal.ii200i : 
Rothberg & Joseph 2006a b; Hopkins et al. 2007d). 
Observations have similarly linked mergers to at least 
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of the quasar p opulation dSanders et alj 119881: 
Canalizo & Stockton 200U iGuvon et all 120061: iDasvra et all 
2007e iBennert et ah i2008l) . Although the precise role of 
mergers is debated, the existence of tight correlations 
between black h ole mass and spheroi d properties such 
as st ellar mass (i Magorrian et al.l Il9 98|). velocity d isper- 
sion jFerrarese & Merritt 2000; Gebh ardt et all 120001). and 
bindi ng energy tHopkins et al.. ,2007Qbl : lAUer & Richstond 
120071; lYoungeretal.ll2008al) imply that the gr owth of black 
holes, dominated by bright quasar phases (Soltan '198 
Saluc ci et al. 1999; Yu & Tremaine 2002; Hopkins et aL 
lOOli. [Shankar etal.ll2009h . is fundamentally linked to the 
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growth of spheroids. 

Despite the importance of galaxy mergers, the galaxy- 
galaxy merger rate and its consequences remain the subject of 
considerable theoretical and observational debate. Halo-halo 
merger rates have been increasingly well-determined with im- 
provements to high-resolution dark matter only simulations, 
with different groups an d simulations yielding increasingly 
£onsistent results (see e.g.'Gottlober et al.'l2001^ 'Stewart et alJ 
2008, 2009a; Fakhouri & Ma 2008b; Wetzel et al. 2009a; 
Genel et al.ii2008 ). But mapping halo-halo mergers to galaxy- 



galaxy mergers is non-trivial, and there are a number of ap- 
parent disagreements in the literature (both theoretical and 
observational) over the absolute rate of galaxy mergers as a 
function of galaxy mass and merger mass ratio. Moreover, al- 
though most of the literature has focused on the most violent 
events (major mergers), various high-resolution simulations 
have shown that a sufficiently large number of minor merg- 
ers (in a sufficiently short period of time) can do as much 
to build bulg mass a s a smaller number of major mergers 



jNaab & BurkertI [20031: iBournaud et al."2005l; I Younger et al.1 
|2008a; Cox et aQlool [Hopkins et al. 2009bi So it remains 
a subject of debate whether or not minor mergers are impor- 
tant for (or may even dominate) bulge formation. Indeed, sev- 
eral questions arise. What is the galaxy-galaxy merger rate as 
a function of mass ratio, galaxy mass, and redshift? Which 
mergers - major or minor - are most important to bulge for- 
mation? Is this a function of galaxy mass and/or bulge-to-disk 
ratio? What is the typical merger history through which most 
of the bulge mass (and, by implication, black hole mass) in 
the Universe was assembled? What room does this leave for 
secular or non-merger related processes? 

The relative importance of mergers of different mass ra- 
tios is critical for understanding the structure of spheroids 
and the related processes above. Although, in principle, 
a sufficiently large number of minor mergers could build 
the same absolute bulge mass as a couple of major merg- 
ers, various simulations have shown that the two scenar- 
ios produce very different structural properties in the rem- 
nants, including rotation and higher-order kinematics, flat- 
tening and isophotal shapes, profil e shapes, central densi- 
ties, e ffective radii, and triaxialities (Weil & Hernauist 1994, 
19961: iBurkert et al.. ,2008; ,Naab et al. 2009; Hoffman et al. 



20091) . Sufficiently minor mergers may not even preferen- 
tially form elliptical-like "classical" bulges, but rather disk- 
like "pseudo-bulges," which have typi cally been associated 
with secular (non-merger) processes ( Younger et al.l l2008at 
lEliche-M oral et al. 2008). These different channels clearly 

^ In this paper, we take "bulge" to refer to cla.ssical spheroids (Sersic in- 
dex > 2, somewhat dispersion-supported spheroids), which are believed to 
primarily form in mergers, unless otherwise specified. We also use the term 
"bulge" to refer to any classical spheroid - small bulges in disks through SO 
and elliptical galaxies. 



imply dramatically different formation timescales and histo- 
ries, also important for understanding the star formation histo- 
ries, stellar populations, colors, abundances and a-enrichment 
of spheroids, and their gradients (see e.g. Mihos & Hernauisl 
ll994at lHop"kins et al.l2009ilRuhland et al.l2009tlFoster et all 
12009). 

Basic questions such as whether or not individual galax- 
ies move continuously in small increments in the Hubble di- 
agram, or can change types significantly, depend on the kind 
of mergers in which they form. Clearly, if e.g. small bulges in 
late-type disks are preferentially formed in early major merg- 
ers that destroy the disk (a new disk being later accreted), or if 
they form in situ in minor mergers that only partially affect the 
disk, the implications for their evolution and the demograph- 
ics of bulges and disks are substantial. Moreover, to the extent 
that starburst and/or AGN activity are coupled to bulge forma- 
tion in mergers, the magnitude, duty cycles, and cosmological 
evolution of these events is clearly directly linked to the kinds 
of mergers triggering activity - one would expect continuous, 
high-duty cycle but low-level activity from sufficiently mi- 
nor mergers, with more dramatic, dusty, shorter duty-cycle 
activity in major mergers ( Hop kins & Hernquistll2006l 120091 : 
Hopkins et al. 2009d). Whether or not observations could, 
in principle, see evidence of merger-induced activity in these 
systems also depends on the kinds of mergers that dominate 
bulge formation, as does the question of whether or not ev- 
ery bulge/massive elliptical passed through a ULIRG/quasar 
phase in its formation. 

Conditions are ripe to address these questions. In addition 
to the convergence in theoretical predictions of the halo-halo 
merger rate, observations have greatly improved constraints 
on halo occupation statistics: namely, the stellar mass distri- 
butions of galaxies hosted by a halo/subhalo of a given mass. 
Observational constraints from various methods yield consis- 
tent results with remarkably small scatter (discussed below), 
and have been applied from redshifts z = — 4 (albeit with 
increasing uncertainties at higher redshifts). Furthermore, to 
lowest order, many of the most salient properties for this ap- 
plication appear to depend primarily on halo mass. That is, 
at fixed Mhaio other properties such as redshift, environment, 
color, satellite/central galaxy status, or morphology have a 
minor impac t (see e.g. Yan et alJ[2003t IZentner et al J 12005 : 
Tinker et al J 120051: ICooravl I2005t iWeinmann et all I2006at 
van de n Bosch et al.1 120071; Iconrov et al.1 120071: IZheng et alj 
2007), the dependence of these properties on mass and one 
another all being expected consequences of a simple halo oc- 
cupation distribution. Populating well-determined halo merg- 
ers with well-constrained galaxy properties can yield predic- 
tions for the galaxy-galaxy merger rate without reference to 
any (still uncertain) models of galaxy formation and wit h 
small uncertainties (as demonstrated in Stewart et alj|2009^ . 
Meanwhile, numerical simulations are beginning to converge 
in predicting how the efficiency of bulge formation scales with 
merger mass ratio and other basic parameters (orbital parame- 
ters, gas fractions, etc.), making it possible to robustly predict 
how much bulge should be formed by each event in a given 
merger history. 

In this paper, we present a simple, empirical model using 
this approach to predict galaxy-galaxy merger rates, and the 
relative contribution of mergers as a function of mass ratio, 
each as a function of galaxy mass, redshift, and other proper- 
ties. We compare variations to the modeling within the range 
permitted by theory and observations, and show that the ques- 
tions above can be answered in a robust, largely empirical 
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fashion, and without reference to specific models of galaxy 
formation. 

In §|2]& §[3] we outline the semi-empirical model adopted 
and show how observed halo occupation constraints lead to 
galaxy-galaxy merger rates. In §|4]we use this semi-empirical 
model to predict galaxy-galaxy merger rates as a function of 
mass ratio, galaxy stellar mass, and redshift. In § |5] we use 
these predicted merger rates together with the results of high- 
resolution simulations to identify the relative importance of 
different mergers as a function of mass ratio, mass, and red- 
shift. In §|6] we extensively vary the model parameters to test 
the robustness of these conclusions. In §|7] we summarize our 
conclusions. 

Except where otherwise specified, we adopt 
a WMAP5 cosmology with (17m , , a s, «.,)= 
(0.274, 0.726, 0.705, 0.812, 0.96) (KomatsVeLal] 120091) 
and a Chabrier (2003) IMF, and appropriately normalize 
all observations and models shown. The choice of IMF 
systematically shifts the normalization of stellar masses 
herein, but does not otherwise change our comparisons. 

Throughout, we use the notation Mgai to denote the bary- 
onic (stellar-Hcold gas) mass of galaxies; the stellar, cold gas, 
and dark matter halo masses are denoted M,, Mgas, and Mhaio, 
respectively. When we refer to merger mass ratios, we use the 
same subscripts to denote the relevant masses used to define 
a mass ratio (e.g. /igai = Mgai.i/A^gai. i), always defined such 

that 0<fl<l (Mgal, 1 > Mgal,2). 

2. THE SEMI-EMPIRICAL MODEL 

In order to track merger histories with as few assumptions 
as possible, we construct the following semi-empirical model, 
motivated by the halo occupation framework. Essentially, 
we assume galaxies obey observational constraints on disk 
masses and gas fractions, and then predict the properties of 
merger remnan ts. The model is described in greater detail in 
iHopkins et al.l (i2009g), and a simi lar variant bas e d on sub - 
halo mergers presented in detail in IHopkins et al.l (l2008dlbh . 
but we summarize the key properties here. 

Note that what we describe here is our "default" model. In 
§ |6] we will systematically vary every part of this model, and 
show that our conclusions are robust. 

2.1. Step 1: The Halo+Subhalo Mass Function 

At a given redshift we initialize the haloH-subhalo popula- 
tion, onto which we will "pa i nt" ga laxies. We do this follow- 
ing the standard [Sheth et aP (1200 Ih mass function, calibrated 
to match the output of high-resolution N-body simulations. 
Note that the total halo mass function (which is what we use 
to map central galaxy stellar mass to primary halo mass) is 
nearly identical to the haloH-subhalo mass function - as such, 
including the subhalos explicitly at this stage (which we can 
do following the fits or methodology in e.g. Vale & Ostriken 
2006; Gao et al. 2004; Nu rmi et al.l2006tivan den Bosch etall 
2005) makes almost no difference to any of our results. This 
choice, and the known ^5% uncertainties in halo mass func- 
tions are negligible compared to other uncertainties in our 
model. We will also, in § |6] consider how systematic changes 
in the halo MF owing to different cosmological parameters 
affect our results, and show that the differences are minimal. 

This defines a primary halo mass function. There are of 
course subhalos in each halo - defining the systems that will 
ultimately merge and therefore be of interest here - but these 
come in the first place from mergers of other halos, which we 
will need to model. We describe this below. 



2.2. Step 2: Painting Central Galaxies Onto Halos 

At a given redshift, we use the halo occupation formalism to 
construct a mock sample of galaxies, each with a parent halo 
or subhalo. Specifically, we begin with the observed galaxy 
stellar mass function (MF), which we take as given. Gener- 
ally, we will be interested in all galaxies (i.e. the total galaxy 
stellar mass function). But if we wish to specifically identify 
only gas-rich or star-forming galaxies, we can just use obser- 
vations that separate the mass function of those galaxies alone 
(i.e. that of star-forming or "blue" galaxies)Q 

We then assign each galaxy to a halo or subhalo in a sim- 
ple manner followin g the standard halo occupation method- 
ology described in iConroy & Wechslerl ([2002); ensuring, by 
construction, that the galaxy mass function and galaxy clus- 
tering (as a function of stellar mass, galaxy color, and physi- 
cal scale) is exactly reproduced. This essentially amounts to a 
simple "rank ordering" procedure, whereby both galaxies and 
halosH-subhalos are rank-ordered in maso, and then assigned 
to one another in a one-to-one fashion. 

2.2.1. Uncertainties in Galaxy Masses and the Mapping 

The observational uncertainties in the galaxy abundance are 
a potential source of uncertainty in the model, especially at 
high redshifts. We will therefore consider three different de- 
terminations of the total galaxy stellar MF as a function of red- 
shift, in our "default" model (the model is otherwise identical, 
but will use one or another determination of the MF). Further 
va riations are explored in § |6] as well. These are taken from 
(D iConroy & Wechsled (12009) (our "default" choice, when no 
other is specified), (2 ) Perez- Gonzalez et aP (l2008bl) . and (3) 
iFonta na etall (120041) . 

Each of these stellar MFs is directly constrained as a 
function of redshift out to at least z > 4, more than suffi- 
cient for our (generally lower-redshift) predictions here. The 
specific choices are not important - we could just as well 
use a number of other observations in the literature (e.g. 
I&undv et al. 2005; Pannella et al. 2006; Franceschini et aU 
%06;' Borch et al. 2006; Fontana et al. 2006; Brow n et aTl 
2007; Marchesini et al. 2009; Kajisawa et al. 2009). We 
choose these three because they bracket the extremes in 
all of these different observations. Specifically, the differ- 
ences between mass functions (l)-(3) represent the range 
across all these different samples. Those, in turn, re- 
flect uncertainties owing to a combination of cosmic vari- 
ance, selection/completeness, and different methods used 
to determine galaxy stellar masses (which can easily con- 
tribute ^ 0.3 — 0.5 dex uncertainties in especially the high- 
redsh ift stellar m ass measure ments; see iMoster et al.l l2009t 
!Marchesini et al.l 2009; Behro ozi et al.1 120101) . Thecompi- 

' At z < 2, we adopt the measurements of the type-separated galaxy stel- 
lar mass functions from lllbert et alj 12009^. if such separation is desired. At 
redshifts z > 2, type-separated MFs are no longer available, so we simply 
assume all systems are star-forming; however, the fraction of massive galax- 
ies that are "quenched" and red has become sufficiently low by z = 2 (and is 
rapidly falling) that it makes little difference (e.g. adopting the upper limit - 
that the red fraction at all masses at z > 2 is equal to that at z = 2 - makes no 
difference to our predictions). 

In what follows, the term "mass" of a subhalo will always, unless oth- 
erwise specified, refer to the "infall" or maximum mass of the subhalo, pre- 
accretion (i.e. the maximum mass the system had while it was still a primary 
halo, before being accreted and becoming a subhalo). This is what is impor- 
tant for the galaxy properties, and what enters into the HOD rank-ordering 
methodology. Obviously, at a later time after accretion, tidal stripping can 
arbitrarily reduce the subhalo mass, but this should not change the original 
accreted galaxy mass. 
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lation in iPerez-Gonzalez et alj (I2008H) draws from a large 
number of different observations, and th us is a representative 
"average." The mass function choice in IConrov & Wechsled 
(|2b09) is interesting because the authors do not directly 
present a stellar MF measurement. Rather, they present a fit- 
ted stellar MF versus redshift that is fitted to both various ob- 
servations of stellar masses and to the history of galaxy SFR 
versus mass and redshift. As such, the well-known discrep- 
ancy between high-redshift SFUs and s tellar mass buildup is 
apparent (see 'H opkins & Beacoml20 06l). The authors address 
this by adjusting the stellar MF at high redshifts to fit the SFR 
measurements. As such, it represents a completely indepen- 
dent constraint on the stellar MF evolution, and is representa- 
tive of the stellar MF that would be obtained with fairly rad- 
ically different assumptions (for example, a stellar IMF that 
evolves with redshift). To this extent, it provides something 
of an upper limit to the uncertainties in the galaxy mass-halo 
mass mapping. 

Note that other methods of HOD fitting (other than our 
rank-ordering approach) yield very similar results. Allow- 
ing e.g. for scatter in the Mgai — Mhaio relation or fitting some 
prior assumed functional form (say e.g. a double-power law 
relation between these quantities, separately for central and 
satellite galaxies, with assumed lognormal scatter), yields 
little difference from the rank-ordering method an d very 
small (< 0. 1 dex) scatter (see e.g. [C onrov & Wechsl ed 120091: 
Wang et al. 2006; Perez- Gonzalez et al.i i2008aj: iBrown et al.l 
2008; Behroozi et al. 201(3). We consider such methodologi- 
cal distinctions in §|6] and find that they make little difference 
to our conclusions. 

2.2.2. Attaching Other Galaxy Properties 

If and when other galaxy properties are required as in- 
put for the model, these too a re as signed according to ob- 
servations. For example, in § 15.41 we consider galaxy gas 
fractions. These are assigned to each galaxy according the 
observed correlations between galaxy gas mass and stel- 
lar mass, which have been quantified at a range of red- 
shifts from z = — 3. As discussed there, we have com- 
piled observations from the available sources, spanning this 
redshift range and a stellar mass range from ^ 10'° — 
10'^ Mq (more than sufficient dynamic ran ge for the pre- 
dictions of interest here), specifically from 'Bell & de Jons' 
(12000); McGaugh (2005); Calura et al. (2008); Shaplev et al. 
(2005); Erb et al. (2006); Puech et al. (2008); Mannucci et al. 
J2009); Crescietal. (2009); Forster Schreiber et al. (2009); 
|Erb (2008); Mannucci et al. (2009). 

We find that these observations can be well-approximated 
as a function of galaxy stel lar mass and redshift with the fit- 
ting functions presented in iStewart et al.l (I2009bl) . and there- 
fore simply adopt those, with a simple lognormal scatter 
term of ^ 0.2 dex representative of the observed scatter at a 
given stellar mass. But adopting any individual measurement 
of these quantities instead gives very similar results^ Of 

" At z = 0, t he gas fractions meas ured are based on measured atomic 
HI gas fractions; 'Bell & de Jons' (WoV) correct this to include both He and 
molecular Ht; McGaugh (2005) correct for He but not H2; Kannappan 
(3)04) gives just the atomic HI gas fractions (this leads to slightly lower 
estimates, but still within the range of uncertainty plotted; H2 may account 
for ~ 20 — 30% of the dynamical mass, per the measurements in Jogee et aJJ 
120051) . We emphasize that these gas fractions are lower limits (based on ob- 
served HI flux in some annulus). At z = 2, direct measurements ai'e not 
available for most samples; the gas masses from Erb et al. (2006) are esti- 
mated indirectly based on the observed surface densities of star formation 



course, the uncertainties in gas fractions and other properties 
increase at higher redshifts, but these are still sub-dominant in 
their ultimate effects compared to the growth in uncertainty in 
the stellar mass function itself. 

2.3. Step 3: Halo-Halo Mergers 

The next step towards identifying merging galaxies is to 
identify merging halos. The halo-halo merger rate (i.e. the 
rate at which formerly primary halos are accreted onto other 
primary halos and thus become subhalos) is well-measured in 
cosmological simulations and defined in the extended Press- 
Schechter formalism. In our default model, we determine 
this rate, as a function of halo mass, merging halo mass ra- 
tio, and redshift, from the fits presented in Fakhouri & m3 
(2008l j), determined from high-resolution N-body simula- 
tions (Springel et al.P2005c). We can include scatter in this as 
well; Fakhouri & Ma (2008a) quantify the scatter in halo-halo 
merger rates across populations. For a given halo population 
and arbitrary time interval, we can then statistically assign all 
halo-halo mergers that have occurred in that interval. We vary 
the determination of the halo-halo merger rate in § |6] (using 
e.g. the results of different dark matter simulations, simula- 
tions with gas included, and selecting halos differently), and 
find that it makes little difference. 

The galaxy properties of the secondary are defined at the 
time when it first becomes a subhalo - i.e. when the secondary 
halo mass is its infall/maximum pre-accretion mass (when the 
halo-halo merger occurs). The primary galaxy continues to 
obey the normal Mgai — Mhaio relation for its total halo mass. 
So if the Mgai ^ A^haio relation evolves with redshift, the pri- 
mary galaxy will lie on the relation defined at the moment the 
actual merger occurs. The secondary will lie on the relation 
defined at the moment of the halo-halo merger (the subhalo in- 
stantaneous mass by the time of merger being much smaller). 
This is the basis for the rank-ordering method of assigning 
subhalo populations, and is key to the agreement between the 
observed small-scale clustering of galaxies and that produced 
by the models here. 

2.4. Step 4: From Halo-Halo to Galaxy-Galaxy Merger 

Of course, a halo-halo merger is not a galaxy-galaxy 
merger We need to follow the recent subhalos until the time 
when the galaxies themselves will actually merge. There are 
two general methods to do this. The first, and our choice 
in this "default" model, is to assign a "merger time" - es- 
sentially, a delay timescale, usually calibrated from the re- 
sults of high-resolution simulations, that represents the time 
for orbital decay from the initial halo-halo merger. This is 
the standard approach adopted in many semi-analytic models 
and models based on the extended Press-Schechter formal- 
ism. There are different possible choices for this merger time, 
and we will consider several of them in § |6] For our default 
model, we adopt the most common: the dynamical friction 
time. Specifically, we use the dynamical friction timescale 
for galaxy-galaxy merger with respect to the initial halo-halo 

and assuming that the z = Kennicutt law holds. However other observa- 
tions suggest that it does indeed hold iBouche et al. 2007). and other indi- 
rect estimates yield similar results dCresci et al. 2009; Mannucci et al. 200'|). 
Moreover recent observations have been able to directly measure the molec- 
ular gas content of galaxies at z ~ 2 — 4, and where available these measure- 
ments agree very well with the assumed /gas (M» , z) relations used here (see 
e.g. Ta cconi et alj.2010.) . We therefore conclude that although appropriate 
caution is due at high redshift. radical departures from our assumptions are 
unlikely. 
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merger, calibrated as a function of galaxy mass, mass ra- 
tio, redshift, ar id orbital parameters in hig h-resolution merger 
simulations in lBovlan-Kolchin et alj (120081) . Specifically the 
formula they present is: 



fdf = 0.0216//" 



(Mi/Ma)'-^ 
ln(l +M1/M2) 



exp[1.977] 



(1) 



where H is the Hubble constant at z, Mi and M2 the halo 
masses at the moment of the halo-halo merger, 77 = // (E) is 
the standard "circularity parameter" {(rj) « 0.5), and rc{E) 
is the circular radius for an orbit with energy E (trivially 
related to the pericentric passage distance rperi//"vii-, more 
usually quoted). Note that this explicitly depends on the 
merger orbital parameters. This allows us to incorporate 
the scatter in merger times seen in full cosmological sim- 
ulations. Specifically, the cosmological distributio n of the 
orbital paramet ers 77 and rpeii / rvir are presented in iBensonI 
(|2005); Khochf ar & Burke rt (2006). Drawing randomly from 
these distributions, we can thus determine some Monte Carlo 
merger population^ 

In general, varying the prescription for the merger "delay," 
across the entire physically plausible range, as we do in § |6] 
leads to factor ~ 2 systematic differences in the merger rate. 

An alternative approach to following galaxies to their 
galaxy-galaxy merger from the halo-halo merger is to track 
the subhalos directly in cosmological simulations. In other 
words, given some cosmological population of halos, we can 
follow the bound substructure of the subhalo after an initial 
halo-halo merger, until the subhalo reaches some sufficiently 
small radius or is completely disrupted, at which point we 
define the "merger" to have occurred. Given this methodol- 
ogy, we can define a "subhalo" merger rate - i.e. a merger 
rate of subhalos being destroyed in primary halos (as opposed 
to a rate of those subhalos simply "falling into" those halos). 
If the simulation is sufficiently high resolution, the subhalo 
final merger/destruction should correspond closely to the ac- 
tual galaxy-galaxy merger between the primary halo central 
galaxy and the satellite galaxy hosted by the subhalo (actu- 
ally, so long as the subhalos are tracked long enough that the 
"remaining" merger time, after subhalo destruction, is small 
relative to the Hubble time, this is sufficient). 

With such a measured subhalo destruction rate from sim- 
ulations, we can convolve with the determined Mgai — Mhaio 
distribution and define directly the galaxy-galaxy merger rate. 
Again we stress that Mhaio for the subhalos, as enters into 
this calculation and determines Mgai is the infall or maxi- 
mum pre-accretion mass (as it generally should be, since the 
galaxy, pre-accretion, would lie on the Mgai — A^haio relation 
for a normal central galaxy, not "knowing" it would be ac- 
creted at some time in the future), as the subhalo instanta- 
neous (post-accretion, stripped) mass goes to zero at the time 
of destruction/merger. Such an approach automatically ac- 
counts for e.g. the orbit distribution of subhalos; however, it 
has its own uncertainties - subhalos are not baryonic galax- 
ies, and the subhalo-subhalo merger time can differ from the 
galaxy-galaxy merger time by factors of up to several (owing 
to the resonant effects of baryons, and finite resolution lim- 
its). In any case, in §|6] we will consider several such deter- 

Note that recently. IWetzell (20101) has .shown that these distributions de- 
pend non-trivially on mass and redshift. However, the sense of that depen- 
dence is such that, if anything, the average merger timescales adopted here are 
always upper limits. As such, incorporating this more detailed dependence 
only strengthens our conclusions (see our discussion in §|6j. 



minations of the subhalo merger rates instead of our merger 
"delay" approach , with each using shghtly different method- 
ologies (see e.g. IStewart e t al."2009a^, 'Krav tsov et al.1l2004t 
IZentner et al.l2005[ 1 van de n Bosch et al. 2005j), and show that 
they yield similar results to our merg er "delay" approach. 
We should note that the formula from Boylan-Kolchin et aT] 
(2008) is calibrated as the total time from halo-halo merger 
to galaxy-galaxy merger in live simulations; thus, it implicitly 
includes all the effects seen in a full simulation (e.g. contin- 
uous mass loss/stripping of the satellites, resonant effects on 
near-passage of the galaxies, and baryonic effects on the ha- 
los). 

Thus, after the convolution with the merger timescale, we 
obtain the rate of galaxy-galaxy mergers as a function of 
galaxy mass M», redshift z, and galaxy-galaxy baryonic mass 
ratio /igai. 

2.5. Step 5 (Optional): Linking Populations Across Different 

Epochs 

If we desire only instantaneous quantities (e.g. the merger 
rate), then what we have already described is entirely suffi- 
cient, and we can simply re-initialize the model at any redshift 
where we want to make predictions. 

However, we will occasionally desire integral quantities 
(for example, how many mergers a given galaxy is likely to 
have experienced in its history). There are several choices of 
method to calculate these. 

The simplest is our choice in the "default" model. Recall, 
we only need to integrate quantities such as the merger his- 
tory statistically. And we will quote quantities such as the 
number of mergers for the primary branch of the halo merger 
tree. As such, we can use the fact that for a halo of mass Mq at 
z = 0, the median primary progenitor mass at higher redshift 
{Mi,{Mo, z)) (and its scatter) is a well-measured function from 
cosmological simulations (in other words, for a halo of mass 
Mo, we know the distribution of progenitor masses M/,(Mo, z) 
at some higher z). We can then integrate over average histo- 
ries, for example for the merger rate: 



(A^merger(Mo,Ai)) 



dz 



;(M,,[Mo,z],z,M)dz. (2) 



In our default model, we adopt the fits to the average halo 
growth tracks/progenitor ma ss distribution given a s a function 
of halo mass and redshift in'N eistein et al.l (120061) . calibrated 
to match the results of high-resolution A^-body simulations. 
Of course, at each redshift z in the integral above, our steps 
1—4 are implicit in obtaining the relevant galaxy properties. 

Note that this is effectively the same as beginning with 
some population of halos and evolving them forward along 
the average growth tracks defined by M/, above, where we in- 
tegrate the halo mass forward using the merger rates and as- 
sign any shortfall to "unresolved" or "diffuse" accretion (since 
only mergers with mass ratios /igai > 10 are of real signif- 
icance here, it makes no difference if this mass technically 
comes from very small halos or truly diffuse material). This 
guarantees that the haloH-subhalo population matches the to- 
tal halo mass function at all times. Re-assigning quantities 
such as Mgai (Mhaio) at each time, according to the observa- 
tional constraints at each redshift, is equivalent to assigning 
the galaxy some net growth in stellar mass (star formation 
rate) and gas mass (inflow minus outflow), but without aprior 
on how much of each is in a bulge or disk componentOj We 

If, for example via some rapid mergers, a galaxy exceeds the 
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can therefore compare the impHed growth rate to observa- 
tions, su ch as the observed stellar mass-star formation rate 
relation dNoeske et al.ll2007al) . For example, the interpolation 
between HODs at different redshifts implies that the SFR in 

(0 3 

Star-forming galaxies scales very crudely as cx Ml 

1.5 



mass 



0.6) ^ 

at z < 2 (with a sharp decline above the turnover 
L,), which agrees well with a number of different ob- 



servational esti mates (see e.g. Blain et al. 1999; Noes ke et al 
2007b a; Papovich et al.|r2006; Martin et al. 2007; B ell et al 
20051; iDamen et al. 2009). In fact, a more rigorous compari- 



son shows quite good agreement - this exercise, for a halo oc- 
cupation model that i s effectively equivalent to th e one used 
here, is presented in IConroy & WechsleJ (|2009), who find 
very good agreement both with observations of the integrated 
SFR versus redshift and t he specific SFR in gal axies at dif- 
ferent masses. Also, see iLee et al] (l2009h and IZheng et"an 
(l2007h . who perform a corresponding analysis and obtain sim- 
ilar conclusions. The variation in predictions between our 
different models is comparable to that in the different mod- 
els considered in these papers. (And recall, at least one of 
the mass function fits we adopt throughout is adjusted specif- 
ically to match the SFR versus galaxy mass and redshift rela- 
tion; as such, agreement with these observations is ensured 
by construction). And despite its simplicity, our approach 
effectively guarantees a match to various other halo occupa- 
tion statistics including stellar mass functions and the fraction 
of ac t ive/passive galaxies as a function of mass dVang et al] 
I200I iWeinmann et aHIIOOeallH: iGerke et al.ll2007h . 

2.6. Step 6 (Optional): The Effects of Mergers on 
Morphology 

Finally, we can couple the model to what is seen in high- 
resolution simulations of galaxy-galaxy mergers, to say what 
effects mergers will have on the galaxy morphologies. 

In the model here, the galaxies are initially (at high red- 
shift) disks]_whena merger occurs, we use the model results 
from iHopkins et al.l (I2009bl) to determine how much of the 
galaxy is converted from disk to bulge. The models used a 
suite of several hundred high-resolution hydrodynamic simu- 
lations, including star formation, black hole growth, and feed- 
back from both to quantify the efficiency of bulge formation 
as a function of merger mass ratioQ orbital parameters, gas 

'^gal(^halo) relation, our model implicitly assumes its growth stalls until the 
halo "catches up"; if it falls below the relation, it experiences the necessary 
gas accretion and star formation for itself to catch up. 

'* It is important to use the "appropriate" mass ratio, for which the 
scalings presented in Hopkins et al. ( 2009b) are calibrated. In detail, the 
authors find that the most dynamically relevant mass ratio is not strictly 
the baryonic galaxy-galaxy ratio fj,g^\ nor the halo-halo ratio /^halo- Rather, 
the important quantity is the tightly bound material that survives stripping 
to strongly perturb the primary. Generally speaking, this is the baryonic 
plus tightly bound dark matter mass (the central dark matter mass, being 
tightly bound in the baryonic potential well, is robust to stripping in simula- 
tions; Ouinn & Goodman 1986; Benson et al. 2004; Kazantzidis et al. 2008; 
IPurcell et al.li20 09). This can be reasonably approximated as the baryonic 
mass plus dark matter mass within a small radius of one NFW scale length 
(rj = fvir/c, where c 10 is the halo concentration; i.e. a few disk scale 
lengths). We find that around this range in radii, our results are not very sen- 
sitive to the precise definition, and in general, the baryonic mass ratio /ig^i is 
a good proxy for the mass ratio calculated with this baryonic+tightly bound 
dark matter mass. However, in e.g. dark-matter dominated systems such as 
low-mass disks, the difference is important in particular for the absolute effi- 
ciency of bulge formation. Therefore, since this is what the simulation results 
are calibrated for, we adopt this mass ratio definition to calculate the dynam- 
ics in a given merger. However, this is not observable; we therefore present 
the predicted merger rates and their consequences in terms of an observable 
and easily-interpreted ratio ^gal (this also makes it possible to compare to 



fraction, and other properties. The large suite of simulations 
spans the parameter space of interest, so there are good sim- 
ulation analogues to the cosmologically anticipated mergers 
here, for which we can simply adopt the fits therein to model 
their bulge formation. 

The details are presented in lHopkins et al.l (l 200 9b|) - specif- 
ically they provide fitting functions (their Equations 7-10 & 
27) that give the exact fraction of both the stellar and gaseous 
disks converted to bulge, as a function of merger mass ra- 
tio, gas fraction, and orbital pa rameters (all quantities that 
are determined in our model)n3 These scalings are all tested 
against the suite of simulations therein, and shown to pro- 
vide accurate fits with < 0.3 dex scatter over the entire dy- 
namic range of merger types of interest for this paper (see 
their Figure 16). In detail, the authors there define "bulge" 
as the stellar population that is dynamically supported by dis- 
persion (as compared to "disk" stars supported by rotation, 
and all cold gas, hence star formation, which is part of the 
disk). They show that this agrees fairly well with decom- 
position of the observable stellar mass profile into an expo- 
nential disk and r'/^-law bulge. But the bulges thus defined 
are generally "classical" bulges, by various observational def- 
initions; as such, this is the "bulge mass" that we predict in 
this paper (pseudobulges, in any case, are believed to form 
primarily via disk instabilities, which we are not modeling 
here). We refer to that paper for the full equations, but in 
an approximate sense, in a merger of mass ratio /i, a frac- 
tion ^ fi of the primary stellar disk is violently relaxed and 
adds its low-density material to the bulge, and a fraction 
~ ( 1 — /gas) A* of the gas loses its angular momentum and par- 
ticipates in a nuclear starburst, adding high-density starburst 
material to the bulge. The factor (1 — /gas) represents the fact 
that the gas cannot be violently relaxed as stars are (because 
it is collisional); rather, it must lose angular momentum to 
the stars (which have a mass fraction in the disk oc 1 — /gas) 
- so in systems with increasing gas fractions, the efficiency 
of gas angular momentum loss decreases. This has very im- 
portant consequences both for understanding what has now 
been seen in essentially all simulations of sufficiently high 
gas-fraction disks (both high-resolution and cosmological; 
see e.g. [S pringel & Hernauist 2005; Robertson et al. 200ot 
Cox et a ll^OOS; Robertson & Bullock 2008; Okamoto et aU 
2005H Scannapieco et al. 2008; Governato et al.l l2007ll20()8lK 
and for the global census 



lus of bulge and disk popula- 
thin disks ( Stewart et al] l2009bl; 



tions and survival of e.g. 
Hopkins et al. 2009g). The true scalings are more detailed 
(see H opkins et al. 2009b) but this represents the important 
physics. 

In § |6] we discuss the consequences of using other, less 
accurate approximations to the behavior seen in full high- 
resolution simulations of galaxy mergers. We find that so long 
as the key scalings with mass ratio are at least qualitatively 
similar, our conclusions are unchanged. 

2.7. Summary 



other results in the literature). Again, the qualitative scalings are the same, 
but it is important to use the full information available in the model to cal- 
culate the merger dynamics. We include all mergers above a minimum mass 
ratio ~ 0.0 1, although our results are not sensitive to this limit so long 
as it is small. 

The efficiency of bulge formation does depend on merger orbital pa- 
rameters - namely the relative inclination angles of the merging disks — so 
we simply draw them at random assuming an isotropic distribution of incli- 
nations (allowing some moderate inclination bias makes no difference). 
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Together, these simple assumptions are sufficient to define 
a "background" galaxy population. There are degeneracies 
in the model - however, we are not claiming that this is 
unique nor that it contains any physics other than the num- 
ber and effects of mergers. For our purposes, the precise con- 
struction of the empirical model is not important - our re- 
sults are unchanged so long as the same galaxy mass-halo 
mass relation and gas fraction distributions are reproduced 
as a function of galaxy mass and redshift. The importance 
of all the above is the following. (1) Observed galaxy stel- 
lar mass functions, the galaxy stellar mass and halo or sub- 
halo mass relations, and the distributions of galaxy gas frac- 
tions as a function of stellar mass and redshift are all repro- 
duced exactly as observed, at all redshifts z = — 4 where 
observational constraints exist, by construction in the model. 
(2) Observed galaxy-galaxy clustering, both on large scales 
(the "two-halo" term, where it is primarily a function of halo 
mass) and small scales (the "one-halo" term, where it reflects 
subhalo/satellite halo occupation statistics), at those redshifts 
z = — 4, are also reproduced by construction. Since this 
is built into the model exphcitly, we do not show such a 
comparison, but'Con roy et all (^2006^ ■ IWang et"aD (^2006). and 
IZheng et al. (2007) all show illustrations demonstrating that 
the simple methodology here yields (again by construction) 
excellent agreement with observed galaxy-galaxy correlation 
functions from scales of 50kpc through lOMpc, at redshifts 
z = 0, 1,2,3 and 4, as a function of galaxy mass or luminos- 
ity, and galaxy color. (3) The only purpose of simulations 
ultimately, in this model, is (as in all halo occupation mod- 
els) to provide a means of linking (statistically) galaxies at 
two different times. In other words, knowing the observed 
distribution of galaxy stellar masses and their separations, we 
use the dynamics which can be followed in the simulation to 
say which of these galaxies will merge in a given time inter- 
val (or again, some statistically approximation such as which 
fraction of the systems inside some distance). Likewise, we 
can ask where the systems that have such mergers "end up" in 
halo or stellar mass. 

Together, this defines the suite of models adopted here. This 
is, by construction, a minimal model, and may leave out im- 
portant details. However, in § |6] we systematically vary the 
model assumptions, and find that our conclusions are robust. 

3. THE CONSEQUENCES: HOW HALO OCCUPATION STATISTICS 
CHANGE GALAXY-GALAXY MERGER RATES 

Figure [T] illustrates how the combination of halo-halo 
merger rates and the observed halo occupation distribution de- 
termines galaxy-galaxy merger rates. 

First, we show the halo occupation function itself (top left 
of Figure [T]): for our purposes, this function is summarized 
in the most important quantity, the average galaxy baryonic 
mass hosted by a halo of a given mass Mgai(Mhaio)- If galaxy 
formation were efficient this would simply be Mgai — /^Mhaio, 
where fb is the Universal baryon fraction, and galaxy-galaxy 
mergers would directly reflect halo-halo mergers. However, 
the relation between Mgai and Mhaio is non-trivial. 

We show several observational constraints on this quantity 
at redshift z = (although we note that it is not redshift- 
independent, according to the observational constraints used 
for the HOD herein). First, the combination of the observed 
abundance and clustering of galaxies of a given mass have 
long been known to set tight constraints on Mhaio(A/gai)- We 
show a recent determination of these constraints, fro m clus- 
tering and abundance of local SDSS galaxies in .Wang et al.l 



( 120061) n Second, we show the empirical "monotonic" or 
"rank ordering" results: it has been shown that good fits to 
halo occupation statistics (group counts, correlation functions 
as a function of galaxy mass and redshift, etc.) over a range 
of redshifts z = — 4 are obtained by simply rank-ordering 
all galaxies and halosH-subhalos in a given volume and as- 
signing one to another in a monotonic one-to-one manner 
(see e.g. IConrov et al.] I2006t IVale & Ostrike^ 120061) . Here 
we plot the results of this exe r cise u sing the z = stellar 
mass functions from 'Bel l et al.l (120031) . Third, we compare 
this to independent esti mates of the averag e Mhaio(A^gai) from 
weak lensing studies in iMandelbaum etalJ (12006). Other in- 
dependent constraints give nearly identical results; these in- 
clude halo mass estimate s in low-mass systeni s from rota- 
tion curve fitting (see e.g. Persic & Salucci 1988;'Persi cet alj 
1996; Borriello & Salucci 2001; Avila-Reese et al. 200l), or 
in high-mass systems from X-ray gas or group kinemat- 
ics (Eke et al. 2004: Yang et aL .2005; Brough et al.. .20061; 
van den Bo sch et al.ll2007h 

Because it is the total baryonic mass, not just the stellar 
mass, that matters for e.g. defining the dynamics in a merger 
(a galaxy can be very massive and, in principle, contain lit- 
tle stellar mass), we will focus throughout this paper on that 
quantity (Mgai and /igai). However, Figure [T] is qualitatively 
identical if we use just the stellar mass instead. In fu- 
ture work (in preparation), we will compare the consequences 
for the HOD, merger rates, and observable quantities such as 
the merger fraction resulting from different choices M*, Mgai, 
Mhaio etc. in the definitions of merger mass ratio. 

The second ingredient in predicting merger rates is the halo- 
halo merger rate determined from A^-body simulations. De- 
fined as the number of halo-halo mergers per primary halo, 
per logarithmic interval in mass ratio /ihaio = Mhaio. 2 /Mhaio i, 
per unit redshift (or per Hubble time), the halo merger rate 
function can be approximated as (,Fakhouri & Ma 2008b) 



Halo dlog/Xhaio dz 



= F(Mha,o)G(z)Mh-a!oexp{(^)'''}. 



(3) 

In these units, halo merger rates are nearly mass and redshift- 
independent; F(Mhaio) ~ 0.03 (Mhaio/ 1 -2 x IOI^Mq)^ ^^ and 
G(z) « (dtJf/dz)"'^^ are weak functions of Mhaio andz, respec- 
tively (in terms of mergers per unit time, this rate increases 
roughly as ( 1 + z )^). A number of other authors gi ve alterna- 
tive fits (see e.g. lStewart et al. 2009a; Genel et alJl200 8). but 
the salient features are similar; weak redshift and halo mass 
dependence (in these units), power law-like behavior at low 
mass ratios with a slope of roughly ji^J^^ and an excess above 
this power-law extrapolation at high /ihaio- 

The (fractional) contribution to halo growth from each in- 
terval is just /ihaio times the merger rate; we plot this quan- 
tity in Figure [T] (top right) as we are ultimately interested in 
which mergers contribute to bulge growth. Because halo-halo 
merger rates go roughly as <iA^,neiger/'^logAihaio oc /i^j^ (re- 
flecting the shape of the halo mass function itself, and the 
nearly scale-free nature of CDM cosmologies), similar mass 
is contributed to the halo from each logarithmic interval in 

Athalo- 

Convolving the halo-halo merger rates with the HOD (i.e. 

These authors and others determine constraints in terms of galaxy stellar 
mass M» ; where necessary, we use our standard fit to observed gas fractions 
as a function of galaxy stellar mass and redshift (see § |2) to convert freely 
between baryonic (Mgai) ™d stellar (M*) galaxy masses at all redshifts. 
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1 . — Top Left: Halo occupation: median baryonic galaxy mass of the primary galaxy in a halo or subhalo, as a function of that halo or subhalo's mass (at 



z = 0). For central galaxies, this refers to the total (primary) halo mass; for satellite galaxies, to their subhalo "infall" mass (i.e. mass the now-subhalo had at the 
last time it was a primary halo). D otted line represe nts maximally efficient s tar formation (/;, is the Universal baryon f raction). We compar e empirical constraints 
from clustering W ang et"ani200 6). weak lensing I Mandelbau m et ani200g ). and abundance matching Cmonotonic MConrov et al.ll200(il) . Our default model is 
constructed to match these constraints. Top Right: Differential contribution to growth from different mass ratio mergers, i.e. merger rate per logaiithmic interval 
in mass ratio and unit redshift, dMndgcisdlog/x^' dz^', weighted by ii. Dotted line is for halos or, equivalently, would be for galaxies if Mgui oc Mhaio were 
the actual HOD. Blue solid line is the result for ^ L, (M* fa lO" Mq) galaxies, given the observed HOD. Bottom Left: Cumulative contribution of different 
mass ratio mergers to the z = spheroid mass density (or halo mass density), integrated over all galaxy (halo) masses. Bottom Right: Differential version of the 
same. Because galaxy mass is not simply proportional to halo mass, bulge growth is dominated by major mergers while halo growth is contributed to by a wide 
range of mass ratios. We focus on Mgui rather than just the stellar mass M» because the former, not the latter, matters for the dynamics in mergers; however the 
conclusions using M* are qualitatively identical. 



populating each halo with a galaxy of the appropriate mass), 
and the appropriate time lag between halo-halo and galaxy- 
galaxy merger, we obtain the galaxy-galaxy merger rate, now 
in terms of the galaxy-galaxy mass ratio /i^ai = A^gai.2/A^gai, i- 
In Figure [T] we compare this function (evaluated at ^ L^, or 
Mgai ~ 10^' Mq, where most of the stellar mass in the Uni- 
verse is concentrated) to that obtained for halos. 

Integrating over the galaxy history in our models. Figure [U 
(bottom left) shows the fraction of the total z = bulge/halo 
mass contributed by mergers with a mass ratio above some 
Mgah i-e. 



/bulge (> Aigal) = 



1 



Mbulee 



8(^^„l//.tgal)dmbul; 



(4) 



where refers to the mass ratio of the merger that formed 
each differential unit dmbuige of the final bulge mass Mbuige, 



and 8(x) = 1 for x> \ {p-'^^\ > /igai) and 8(x) =0 for x < 1 
(n'^^i < /igai)- Note that this can be defined over the bulge 
mass of an individual galaxy, over all bulge mass in galaxies 
in a narrow interval in mass Mgai, or over all galaxies (i.e. in- 
tegrating over the bulge mass function). We show the latterFl 
We also show the differential version (bottom right): the frac- 
tion of z = bulge mass contributed per logarithmic interval 
in merger mass ratio d/buige/dlog/igai. 

For halos (or equivalently for galaxies if the trivial map- 
ping Mgai A^haio wcrc truc), the distribution of bulge mass 
fraction from mergers of different mass ratio is quite broad, 
as expected: only ^ 50% of halo mass comes from mergers 

" In principle, some bulge mass could come from redshifts before our 
"initial" tracking of each halo, but in practice at any redshift, most of the mass 
has assembled relatively recently, so our results do not depend sensitively on 
the initial conditions. 
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with /i > 0.1. Because halo-halo merger rates are nearly self- 
similar, the differential version of this reflects the instanta- 
neous rate also shown, with similar contributions per logarith- 
mic interval in halo mass ratio. It is still the case, though, that 
ten 1:10 mergers are less common than a 1 : 1 merger, meaning 
that major mergers dominate (Stewart et al. 2008). 

In contrast, galaxy bulge assembly is biased much more 
towards high-mass ratio mergers, at least for ^ L» systems 
which dominate the mass density. This owes to the nature of 
halo occupation statistics: at low masses, galaxy mass grows 
rapidly with halo mass (galaxy formation is increasingly effi- 
cient as one moves from low masses closer to ~ L,). Upon 
reaching ~ L*, however, star formation shuts down relative 
to halo growth - in terms of the HOD, the scahng of galaxy 
mass with halo mass transitions from steep Mg^\ oc M^^^'^ '^ 
at low masses to shallow Mgai cx at high masses. In 

short, as halos grow in mass past ~ IO'^Mq, galaxy masses 
"pile up" near ~ (Mgai ^ 10" M©), as can be seen in Fig- 
ure [T] Since halo-halo merger rates are nearly self-similar in 
terms of halo-halo mass ratio, the "pileup" of galaxies near 
this mass means that a wide range of halo-halo mass ratios 
will be compressed into a nar row range of galaxy-galaxy mass 
ratios near /tgai ^ 1 (see also lMaller et al.ll2006l: IStewart et al.l 
|2b09a; Stewart 2009). 

This is a general statement: in any scenario where 
A^haio/A^gai has a minimum, the galaxy-galaxy merger rate will 
be weighted more towards major mergers than the halo-halo 
merger rate around (in particular at masses slightly above) that 
minimum. The minimum in Mhaio/A^gai is empirically well- 
established, and occurs near ^ L,, where most of the mass 
density lies. It is therefore inevitable that the contribution to 
the integrated bulge mass will be more weighted towards ma- 
jor mergers than to the halo mass. 

This is also easily understandable in terms of the mass func- 
tions of galaxies and halos. The halo mass function does not 
feature a sharp break, so the mass density of halos is broadly 
distributed over several orders of magnitude in halo mass. In 
contrast, the galaxy mass function reflects inefficient star for- 
mation at low and high masses, with a sharp break, and so 
the galaxy mass density is concentrated in a narrow range (a 
factor 3) around the break L,. The mass of subunits (which 
broadly reflects the global mass function) in halos therefore 
includes contributions from a wide range of mass ratios. In 
contrast, the bulge growth of a galaxy is dominated by sys- 
tems near ~ L*. At masses < a few L,, this means major 
mergers will be most important. It is not until an ^ galaxy 
represents a minor merger (i.e. galaxy masses > 3 L*) that mi- 
nor mergers (again, mergers of those ~ galaxies) begin to 
dominate the mass assembly. 

4. GALAXY-GALAXY MERGER RATES 

4.1. Scaling with Mass and Redshift 

We first examine the galaxy-galaxy merger rate, given these 
empirical constraints. Figure|2]shows the number of mergers, 
as a function of mass ratio, that a typical galaxy of a given 
z = stellar mass has experienced since z = 2. We emphasize 
that this is for a sample mass-selected based on their z — Q 
masses; a sample selected at the same mass at higher red- 
shift will have a systematically larger number of mergers in 
a similar time or Az interval (and will be higher mass by 
z = 0). Below, we discuss how these predictions compare 
with observations; the two generally agree well. We com- 
pare with the number of halo-halo mergers as a function of 
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Fig. 2. — Top: Median number of galaxy-galaxy mergers since z = 2 above 
a given baryonic mass ratio fi„jj for a z = galaxy of the given stellar mass. 
The same mass selection at higher redshift intervals will systematically in- 
crease the number of mergers. The predictions here agree well with observa- 
tional estimates ide Ravel et al. 2009; Conselice et al. 2009; Lin et al. 2008). 
Bottom: Same, for halo-halo mergers. Owing to the shape of the HOD, the 
merger rate is a steep function of galaxy mass around ~ L, , whereas the halo- 
halo merger rate is only weakly mass-dependent. Galaxies near ~ L« have 
had ~ 1 major merger since z = 2. 

halo mass ratio, for typical corresponding halo masses. The 
two are quite different, for the reasons discussed in § [3] es- 
sentially, at low masses, Mgai M^aio' ^° ^ ^-^ Mhaio merger 
becomes a 1:9 /igai merger, and merger rates at each /igai 
are suppressed. At high masses, Mgai oa Mj^^j^,, and rates are 
correspondingly enhanced. At ^ 10" Mq, where most of 
the spheroid mass density of the Universe resides, the typ- 
ical galaxy has experienced ^ 0.5 — 0.7 major (fig^i > 1 /3) 
mergers since z = 2, a fraction that corresponds well to the 
observed fraction of bulge-dominated early-type systems at 
these masses (see e.g. lBell et al.ll2003l) . At most masses (ex- 
cepting the highest masses, where the shape of the HOD 
yields a strong preference towards minor mergers), the to- 
tal number of mergers with /igai ^ 1/10 is a factor ~' 2 — 3 
larger than the number with /igai > 1/3, and at low mass ra- 
tios /igai ^ 1/10, the merger rate asymptotes to a power-law 

•.u AT/ \ -(0.25-0.5) 

With A^(> /igai) OC/igJ 

Figure [3] shows how the median merger rates evolve with 
redshift, for four different intervals in mass. We compare 
the rate of major (/igai > 1 /3) and major-i-minor (/igai > 1 /lO) 
mergers. We also compare different constraints on the HOD, 
from fitting different galaxy mass functions and clustering 
data. Specifically, we show the default model here, where 
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Fig. 3. — Merger rate (number of mergers per galaxy per Gyr) as a function of redshift for different galaxy mass intervals (all Mg^i above some minimum 
baryonic mass in Mq at each redshift). Solid, dashed, and dot-dashed lines correspond to three different halo occupation model constraints: uncertainties are 
small at z < 2. Colors correspond to the range of merger mass ratios. 



the function Mgai(Mhaio) is determined from fits to observed 
clustering, stellar mass functions, and star formation rate dis- 
tributions in IConrov & Wechsler! (|2009); we compare the re- 
sults adopting a monotonic ranking between galaxy and halo 
mass and using the redshif t-d ependent stellar mass function s 
from lFontana etal.1 ( 120061) or iPerez-Gonzalez et all ( l2008bl) . 
Further variations are discussed in detail in § |6] These illus- 
trate the robustness of the model: empirical halo occupation 
constraints are sufficiently tight that they contribute little am- 
biguity in the resulting merger rate at z < 2. Above z — 2, 
the results begin to diverge, as the stellar mass function is 
less well-determined (and few clustering measurements are 
available); however these higher redshifts have relatively little 
impact on the predictions at low-z. The major merger rate in- 
creases with mass with a slope of roughly oc (1 +z)''^~^", but 
this is mass-dependent. The evolution in the galaxy-galaxy 
merger rate is somewhat shallower than the redshift evolution 
of the halo-halo merger rate (which scales with the Hubble 
time) as a consequence of the redshift evolution of the HOD. 

4.2. Comparison with Observations 

In Figure m we compare these predictions to observed ma- 
jor merger fractions. As most measurements of the merger 
fraction do not have a well-defined mass selection, we first 
simply consider a large compilation of observational results 
compared to the predicted merger rate of ~ galaxies. For 
now, because we are considering a range of mass and ob- 
servational methodologies, we simply convert the predicted 
merger rate to an observed merger fraction assuming a con- 
stant observable lifetime fobs, here showing fobs = 0.5 Gyr and 
fobs = 1 Gyr, typical values in the literature. We also compare 



TABLE 1 
Observed Merger Rates 



Reference 



Selection 



Symbol^ 



Kartaltepe et al. (2007) 
Linet al. (2004, 2008) 
Xue t al. (2004) 
De P ropris et al. (2005) 
Bluck et al. (200^ 
Bundv et al. (2009) 
Bell et al. (2006b,a) 



Pairs 



20/!-'kpc 
30/i^'kpc 
20/i-'kpc 
20/!-'kpc 
20/i-'kpc 
20 ft-' kpc 
20/i-'kpc 

Morphology 



blue triangles 
pink circles 
blue circle 
black asterisk 
orange squares 
green stars 
red pentagons 



Conselice et al. (2009) CAS 

Conselice et al. (2008) CAS 

Cassata et al. (20^ CAS 

Jogee et al. (20()8) visual 

Bundv et al. (2005) visual 

Wolfetal. (2005) visual 

Bridge e^al. (2007, 2010) visual 

Lotz etal. 72006. 2008b) Gini-M20 



pink x's 
blue x's 

cyan inverted triangles 
pink triangles 
light green stars 
orange diamonds 
purple squares 
dark green +'s 



' Selection criterion used to identify merger candidates. For pair samples, this 
refers to the pair separation. For morphological samples, to the method used. 
^ Symbol used for each sample in Figuresl4]&[5l 



the number density of mergers versus mass, at a given red- 
shift. The agreement appears reasonable, but there is a large 
scatter in the observations, mostly owing to different selection 
and merger identification criteria. 

We can also compare the predicted integrated number of 
mergers in Figure |2] to various observational estimates. For 
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Fig. 4. — Top: Major merger fraction of L»galaxies. Observations (witli 
the symbol type for eacJi) ai'e listed in Table [TJ for now, we treat them all 
the same. Shaded range corresponds to the predicted merger rate, convolved 
with an observable lifetime tabs = 0.5 — 1.0 Gyr (color denoted the mass ratio 
range, as labeled). Bottom: Integrated merger rate as a function of galaxy 
mass at z = 0.6 (in absolute units, mergers per Mpc^^ per Gyr above some 
minimum mass), w here observ ations span the greatest dyna mic range. Ob- 
servations are from Bell et al.' ('2006a, blue square), Conselice et al.' f200p, 
pink triangles), Bundy et al. (2009, green stars), and Lopez-Sanjuan et al. 
12009b. purple circle). Here, morphological samples are converted to rates 
with tabs = 1 Gyr; pair samples the appropriate timescales for their separation 
(see text). Predicted rates agree with those observed, but there is considerable 
scatter between various observational selection methods. 



example, Ide Ravel et al.l (12009) use pair-selected samples to 
estimate that ~ 20 - 25% of the > lO'^M© population has 
experienced a merger with mass ratio fi > 0.25 since z = 1; 
IConselice et al.l (12009 ^ estimate a similar number of mergers 
for > lO'^M© galaxies since z ~ I and about again as 
many since z = 2, they also find that the number of merg- 
ers increases significantly with galaxy mass; iLin et al.l(l2008l) 
estimate that ~ 54% of (M* ^ 10" Mq) galaxies have ex- 
perienced a /i > 0.25 merger since z— 1.2. All of these pre- 
dictions can be compared with our predictions in Figure |2] 
and they agree well (especially given different methodologies, 
masses, and redshift ranges involved). 

In order to test more strictly, and to take advantage of where 
observable merger timescales have been rigorously calibrated, 
in Figure|5]we restrict our comparison to galaxy-galaxy merg- 
ers identified observationally using a consistent methodology 
and covering a well-defined mass range. 

First, we consider pair fractions: specifically the fraction 
of major (figai > 1/3) pairs with small projected separations 
rp < 20/i^'kpc (often with the additional requirement of a 



line-of-sight velocity separation < 500 km s '), and stellar 
masses - 1 - 3 x lO'^M© or - 0.5 - 2 x 10" Mq. 
For each mass bin, the pair fractions as a function of red- 
shift can be empirically converted to a merger rate using th e 
merger timescales at each radius. iLotz et al.l (l2008al l2009bllah 
specifically calibrate these timescales for the same projected 
separation and velocity selection from a detailed study of a 
large suite of hydrodynamic merger simulations (including a 
range of galaxy masses, orbital parameters, gas fractions and 
star formation rates) using mock images obtained by applying 
realistic radiative transfer models, with the identical observa- 
tional criteria to classify mock observations of the galaxies at 
all times and sightlines during their evolution. For this spe- 
cific pair selection criterion (if we average over the typical 
distribution of mass ratios for mergers selected in this inter- 
val), they find a median merger timescale of fmerger ~ 0.35 Gyr, 
with relatively small scatter and veiy little dependence on 
simulation parameters (±0.15Gyr)nj We use their median 
fmerger to convcrt the obscrvations to a merger rate. Because 
of the weighting over merger mass ratio and orbital parame- 
ters, for which the explicit dependence is presented in these 
papers, we obtain the same result (within the observational 
error bars) if we convolve our predicted merger rates with the 
explicit observable merger timescale as a function of merger 
gas fraction, galaxy mass, redshift, mass rati o, and orbital 
parameters, as presented in iLotz et al.l (l2009bllah . Complete- 
ness corrections are discussed in the various papers; we also 
adopt the standard correction from Patton & Atfield ( 2008), 
calibrated to high-resolution simulations, for the fraction of 
systems on early or non-merging passages (to prevent double- 
counting systems on mult iple passages); bu t this is relatively 
small (20-40%; see also lotz et aAoOS al). 

The advantages of these pair fractions are that: (1) the mass 
ratio can be determined, leading to little contamination from 
minor mergersQ (2) at such small separations, most such 
pairs will eventually merge, and (3) there is little ambiguity 
in the merger timescale, with only a factor ~ 10 — 20% sys- 
tematic uncertainty in the median/average merger timescale 
in high-resolution calibrations (with a 25 — 50% dispersion 
or variation always present about that median, owing to cos- 
mological variation in e.g. the exact orbital parameters). 

Second, we consider morphologically-selected mergers, 
identified on the basis of by-eye classification or automated 
morphological criteria such as the co ncentration-asymmetry 
(CAS) or Gini-M20 planes (see e g. IConselice etall 120031: 
ILotz et al.ir2004l) . iLotz et al.l (l2008al) also attempt to caUbrate 

The merger timescale from simulations at this radius is shorter than 
the time obtained assuming dynamical friction and circular orbits in e.g. an 
isothermal sphere, as has commonly been done (this is assumed in e.g. both 
[Patton et al. 2002; Kitzbichler & White 2008). This owes to two effects: first, 
angular momentum loss at these radii is not dominated by dynamical friction, 
but rather by exchange in strong resonances that act much more efficiently 
(even allowing for e.g. mass loss by the secondary inside the primary halo, 
which is of course included in the high-resolution simulations, this has the 
net effect of significantly accelerating most mergers). Second, by these radii, 
even initially circular orbits have become highly radial, leading to shorter 
merger times. Because of these effects, the remaining merger time at this 
scale depends only weakly on initial conditions or orbital parameters - es- 
sentially, these processes have erased most of the "memory" of the original 
orbital configuration. This emphasizes the importance of using full simula- 
tions with baryonic effects in calibrating these timescales. 

" Note that many older studies adopt the galaxy-galaxy luminosity ratio 
as a proxy for mass ratio. This is not a bad approximation in e.g. numerical 
simulations, but could be subject to bias from e.g. differential enhancement 
in star formation. Obviously it is preferable to use an actual stellar mass ratio 
where possible. Restricting our samples to just studies with stellar masses, 
however, we obtain similar conclusions. 
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Fig. 5. — Top: Predicted merger rate in different stellar mass intervals, compared to that inferred from close pair studies (projected r <20h~ kpc). Observations 
(with point types) are listed in Table[Tl they are converted to merger rates given the calibration of merger times as a function of pair separation from high-resolution 
N-body simulations in Lotz et al. 1 2008a, 2009b a). We show our default model, as well as the "instantaneous" calculation obtained by assuming all galaxies 
lie on the z = (non-evolving) SDSS HOD shown in Figure[T] We also show the results from Stewart et al. l2009a), using different simulation merger rates, 
merger trees, and HOD constraints. Bottom: Morphologically identified merger fraction in the same mass intervals, compared with predictions. The assumed 
timescale for identification as morphologically disturbed is ~ 1 Gyr, again calibrated from simulations in Lotz et al. (2008a). Observations in this case could be 
contaminated by minor mergers; we therefore show both major and major-l-minor merger fractions. The scatter is larger in morphological samples (with some 
probably contamination), but predictions agree within a factor ~ 2. 



the observable timescale for classification of major merg- 
ers via the Gini-M20 criterion, at rest-frame wavelengths 
and masses of the observations. They find an observable 
timescale fobs(Gini — M20) ^ IGyr, and we adopt that here, 
but note that the predicted timescale in this case depends 
much more sensitively on the depth of the observations, the 
waveband adopted, anc l p roperties such as the gas-richness 
of the merging systems^] Moreover, although, by definition, 
this methodology is complete to events that have violently dis- 
turbed the galaxy, the level of disturbance at a given merger 
mass ratio depends on orbital parameters and galaxy gas frac- 
tions, so a fixed level of disturbance does not correspond to 
a fixed merger mass ratio. Some c ontamination from minor 
mergers is likely. I Jo gee et aP ({2009') estimate empirically that 
^ 30 — 40% of their (by-eye) morphologically-identified sam- 
ple represent contamination from 1/10 < /igai < 1/3 minor 



■^^ F or the explicit dependence on these parameters, see ILotz et al] 
I2009bl l3). Note that the timescale of 1 Gyr here is slightly longer than the 
0.6 Gyr estimated directly from observations in Conselice (2009) (although 
within their quoted Icr eiTor bars) and in the original Lotz et al. 1 2008a) for 
radial 1 : 1 mergers. The difference in the latter owes to the dependence of ob- 
servable timescale on gas fraction as calibrated in Lotz et al. ( 2009a) (here we 
adopt a median appropriate for the median gas fractions of the model galax- 
ies), and from the dependence on merger mass ratio as calibrated in Lotz et alj 
('2009b) (where our value here represents a weighting over the mass ratio 
distribution down to mass ratio fi = 1/3 as appropriate for the observations 
here). In any case, the difference is generally smaller than the scatter between 
different observational estimates. 



mergers. We therefore compare the predicted merger frac- 
tion for /igai > 1/3 and /igai > 1/10. Allowing for this range, 
the observations agree well with the predictions. In particu- 
lar, the observations using a calibration of Gini-M20 or CAS 
specifically matched to hi gh-resolution hydrodynamic major 
merge r simulations (e.g. iLotz et all l2008bt IConsehce et al.l 
I2OOI agree reasonably well with the predicted /igai > 1 /3 
fractions. And external, purely empirical indic ators favor sim- 
ilar merger timescales (see e.g. IConselicdl200 9). 

A quantity closely related to the pair fraction on small 
scales is the galaxy-galaxy autocorrelation function (specif- 
ically that on small scales, inside the "one halo term" where it 
reflects galaxies inside the same parent halo). Effectively this 
generalizes the predicted pair fraction from < 20 /1 ^ ' kpc to all 
scales. But recall, the adopted halo occupation-based method- 
ology is designed, by construction, to match the observed cor- 
relation functions as a function of mass. It is therefore guaran- 
teed that the clustering at scales ^ 100 kpc through > lOMpc 
is reproduced as a function of galaxy m ass and redshi ft (for 
explicit illustrations, see e.g. lConroy et al..2006.; Zheng et al.l 
I2007tlmni"ern]l2006h . 

4.3. Analytic Fits 

It is useful to quantify the predicted merger rates with sim- 
ple analytic fitting functions. 

First, consider major mergers. We find that the major 
merger rate (number of /Xgai > 1/3 mergers per galaxy per 
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Fig. 6. — Results of fitting the major merger rate as a function of mass 
and redshift to a function of the form o f Equ ation \5\ Linestyles con'espond 
to different HOD choices, as in Figures 13151 Left: Normalization, i.e. num- 
ber of mergers per galaxy per Gyr at z = 0, for galaxies above a minimum 
mass Mmm- Right: Redshift dependence, i.e. slope /3 given a merger rate per 
unit time oc (1 -j-z)^ , as a function of minimum mass. The mass-dependent 
A{Mt^in) and [H can be approximated with Equations I6I7I with systematic 
uncertainties of 0.3 dex and 0.2, respectively. 



unit time), for galaxies above a given minimum stellar mass 
threshold (M* > M^a), can be well-fitted by the following 
simple function: 



dA^, 



major 



df 



= A(M™n)(l+z 



,/3(M„, 



[per galaxy], (5) 



i.e. a z normalization A (Mmin) and simple power-law scal- 
ing with redshift with slope /3(Mn,in)- Figure |6] shows these 
quantities, fitted to the predictions shown in Figures |3]|5] as 
a function of the mass M^in- The trends discussed above 
are evident: the normalization of merger rates increases with 
mass, and (albeit more weakly), the dependence on redshift 
decreases with mass. This normalization and redshift varia- 
tion can be approximated with the scalings: 



!0.02 



1 



Mo 



Gyr" 



and 



^(M™i„)major « 1.65 -O.lSlog f^^) 

V M{) / 



(6) 



(7) 



where Mq = 2 x 10 Mq is fixed. There is a systematic fac- 
tor ^ 2 uncertainty in the merger rate normalization A (Mmn) 
at all M,nin, considering the range of models discussed in de- 
tail in §|6]below. The uncertainty in /3(Mmin) is illustrated in 
Figure|6] approximately a systematic A/3 ^ 0.15 — 0.20. 

These fits are for major (/igai > 1/3) mergers. To rough 
approximation, the number of mergers as a function of mass 
ratio scales with the approximate form 



CW(> A^gal) 

df 



« A^ea?'^ (1 - Mgai) [per galaxy] 



(8) 



(derived and discussed in more detail in lStewart et al.ll2009al) . 
This is a good approximation as long as the galaxy is within an 
order of magnitude of ^ L*. For the range where this function 
is a good approximation, it implies an approximately constant 
ratio of a factor w 2 (0.3 — 0.4 dex) of major-i-minor (/igai > 
1 /lO) to major (pg^i > 1 /3) mergers. 

More specifically, we can fit a function of the form of Equa- 
tion|5]to the major-i-minor (fig-A > 1 /lO) merger rate of galax- 
ies above some Mj^in, and obtain the best-fit scalings 
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(9) 



/3(M™i„)„ 



1.50-0.25 log 



Mo 



(10) 



with similar systematic uncertainties in both A{M^i„) and 
/3(Mn,in) to the major merger rate. Note that these equations 
should be treated with caution for the most massive systems 
- the simple fitting functions do not extrapolate to arbitrar- 
ily high mass and the direct numerical results (e.g. Figure O 
should be used. 

5. THE RELATIVE CONTRIBUTIONS TO BULGE GROWTH FROM 
DIFFERENT MERGERS 

5.1. Overview 

Figure [7] illustrates how the efficiency of bulge formation 
scales in simulations. We show how the average B/T re- 
sulting from disk-disk mergers scales with mass ratio (ap- 
proximately /igai, but see § |2]l, gas fraction /gas, and merger 
orbital parameters, according to the fits t o the hydrody- 
namic simulations in Ho pkins et al.l (1200 9b'). To lowest or- 
der, as discussed in § [3] the amount of bulge formed (the 
amount of stellar disk of the primary galaxy that is vi- 
olently relaxed, and amount of gas disk that is drained 
of angular momentum and participates in the nuclear star- 
burst) scales linearl\F1 with the mass ratio of the encounter, 
oc /Zgai. This conclusion - ultimately the important state- 
ment for our analysis - has been reached by numerous in- 
dependent simulation studies, adopting different methodolo- 
gies and numerical techniques, and naturally follows from 
the simple gravitational dynamics involved in violent relax- 
ation (see e.g. [ Hernquist 1989; Barnes & Hernquist 1992; 
'Mihos & Hernquist 1994b, 1996; Naab & Burkert 2001 
Bournaud et al. 2005; Younaeretal. 2008a; di Matteo et aQ 
fSoOTtlHopkins et al. 2009bt lCox etalJI2008n . In addition, ob- 
servational constraints on the efficiency of merger- induced 
star formation support these estimates dWoods et alj 120061 ; 



sup p 

iBarton et al.ll2007l; iWoods & GelleiillOOTh . 

At fixed mass ratio, the bulge formed (remnant B/T) can 
vary considerably depending on orbital parameters of the 
merger, in particular the relative inclinations of the disks (pro- 
grade or retrograde). This variation is shown in Figure |7] 
However, in a cosmological ensemble, this will average out. 
Here, we assume random inclinations, but allowing for some 
preferred inclinations amounts to a systematic offset in the 
B/T predicted and will not change our conclusions regard- 
ing the relative importance of different mass ratios for bulge 
formation. 

Another important parameter determining B/T is the 
merger gas fraction. To lowest order, angular momentum loss 
in gas is suppressed by a factor ~ (1 — /gas); as a result, the 
efficiency of bulge formation {B/T expected for a merger of 
a given mass ratio) is suppressed by the same factor. This 
can have dramatic cosmological implications, because /gas is 
a stro ng function of g alaxy m ass; these are discussed in de- 
tail in iHopkins et al.l (12009 d) . Figure [T] shows the observed 
dependence of disk/star-forming galaxy gas fractions on stel- 
lar mass at z = and z = 2; if we assume mock galaxies on 
the z = 2 relation each undergo a merger of a given mass ra- 

^' Note that this refers to the mass fraction which is violently relaxed, thus 
adding to the bulge. Disk heating and resonant processes that cont ribute to 
the thick disk or disk substructure are different. For example, Hop kins et al.l 
i200 8e) show that disk heating in minor mergers is second-order in mass 
ratio;lPurcell et al. (2009) and Kazantzidis et al. (2009) reach similar conclu- 
sions (albeit with slightly different absolute normalization/efficiency). 
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Fig. 7. — Left: Average bulge-to-total ratio B/T (of the classical bulge) resulting from a single merger of mass ratio ~ ^gai between galaxies with gas fraction 
/gas (from simulations in 'Hopkins et al."2009tf). To lowest order B/T scales as cx: ^g^i (1 — /gas)- Error bar shows scatter owing to the cosmological range of 
orbital parameters. Center: Median observed gas fraction and scatter (error bars) for disks of a given stellar mass at z = (blue) from Bell & de Jong l200l|, 
diamonds). .Kannap pan (2004, squai'es), and McGaugh (2005, circles), and at z = 2 (red) from Erb et al. (2006). Solid lines show fits to the median at each 
redshift. Right: CoiTesponding median B/T expected from mergers at z = 2 with primary of a given stellai' mass and mass ratio fj,„^i (given the observed 
(/gas[M« I z = 2])). Suppression of bulge formation by gas-richness is important for the absolute bulge mass formed (especially at low masses), but because it is 
/igai -independent, does not affect the relative contribution of major/minor mergers. 



tio, then the resuhing B/T at each mass is shown. Bulge 
formation will be significantly suppressed by high gas frac- 
tions in low-mass galaxies, giving rise to e.g. a strong mass- 
morphology relation similar to that observed jStewart et al.l 
|2009b; Hopkins et al1 l2009g|) . However, it is clear in Figure|7] 
that the effect of /gas is a systematic offset in B/T, indepen- 
dent of mass ratio. Because in what follows we will generally 
examine the relative importance of mergers of different mass 
ratios (independent of gas fraction), the inclusion or exclusion 
of the effects of gas on merger dynamics makes little differ- 
ence to our conclusions. 

Given these constraints from simulations on the amount of 
bulge formed in a given merger, and the merger rates pre- 
dicted in § m Figure |8] shows the contribution of mergers of 
different mass ratios to the z = stellar mass in bulges (as Fig- 
ure [Til. We show this for galaxies in a narrow range of total 
stellar mass around three different values, and for the entire 
galaxy population (integrated over bulges of all masses). We 
specifically define this as the fraction of bulge mass formed or 
assembled in mergers that were above a given galaxy-galaxy 
mass ratio /igai. Considering different variations to the em- 
pirical and simulation constraints (see § |6]l, about 60 — 70% 
of the globally integrated bulge mass is assembled by major 
mergers /Xgai > 0.3, with another ~ 30% assembled by minor 
mergers 0.1 < /igai < 0.3, and the remaining '-.^ — 10% from 
a wide range of mass ratios /igai ^ 0.01 — 0.1. 

Note that, for massive galaxies (>a few L*), where "dry" 
mergers become an important channel (assembling mass al- 
ready in massive bulges), there is an ambiguity in the fraction 
/buige(> A'gai) "contributed" by mergers above a given /igai. 
Figure [8] illustrates this. We therefore introduce the distinc- 
tion between bulge formation and bulge assembly, terms we 
will use throughout. 

First, formation: we can define /buige(> Mgai) as the frac- 
tion of bulge mass originally /ormeii, i.e. initially converted 
into bulge mass from disk mass (gas or stars) by mergers 
with some mass ratio /igai. In other words, taking EquationH] 
where for each parcel of mass in the final bulge (dmbuige), the 
"contributing" mass ratio /ig^, is defined by /igai of the merger 
that first made the mass into bulge (regardless of whatever 
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Fig. 8. — Top: Cumulative contribution to original /oraia^ion of spheroid 
mass density from mergers of different mass ratios, for galaxies of given stel- 
lar mass at z = (as in Figure[T). "Integrated" curve refers to the integral over 
all bulge masses (net contribution to global spheroid mass density). Bottom: 
Same, but showing the contribution to the integrated spheroid assembly. Ma- 
jor mergers dominate neai' ~ L, ; minor mergers become more important at 
lower/higher masses. Assembly by minor dry mergers (of bulges first formed 
in more major mergers) occurs at the highest masses. 



merger ultimately brings it into the final galaxy). This quan- 
tity answers the question "what kind of merger destroyed the 
progenitor disks of these galaxies?" or "what kind of merger 
created most bulges in the first place?". 
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Fig. 9. — As Figure [8] (assembly), comparing the full result to an instan- 
taneous calculation from convolving the z = HOD of IWang et alj 420061) 
with instantaneous z = merger rates. Because the shapes of the merger rate 
versus ^ihalo and M„i,i(Mi,iiio) functions do not evolve strongly with redshift, 
both yield similar results. 

Second, assembly: we can define /bulge (> Mgai) the frac- 
tion of bulge mass assembled into the main branch of the 
galaxy by mergers with mass ratios > /igai. Here, we take 
Equation |4] where, for each parcel of mass, /ig^, is defined 
by /igai of the merger that brought it into the main progen- 
itor branch of the final galaxy. This answers the question, 
"what kind of merger brought most of the present-day bulge 
together?" or "what kind of merger has affected most of the 
mass in the bulge?". 

Clearly, the two are equal if all bulge mass is formed "in 
situ" in the main progenitor - i.e. if the secondary galaxies 
have no pre-existing bulges. Indeed, at low masses, where 
most galaxies have little bulge, there is little difference. 

At high masses, however, there is a dramatic difference. 
This makes sense: high-mass galaxies grow primarily by 
dry mergers. Galaxies first become bulge-dominated around 
^ L*, then assemble hierarchically. Since most bulge mass 
is first formed around ~ L*, where we see bulge formation is 
dominated by major mergers, we obtain the result that bulges 
in high-mass galaxies are pnmaiily formed in major mergers. 
However, as they grow in mass via dry mergers, minor merg- 
ers become increasingly important to the assembly of the most 
massive systems (minor mergers bring together bulges that 
have already been formed). Both definitions are clearly im- 
portant and have their applications; however, because of e.g. 
the importance for remnant kinematics and growth histories, 
we will generally adopt the latter (assembly-based) definition 
in what follows. 

The most important determinants of Figure[8]are the shapes 
(logarithmic slopes), not normalizations, of the halo merger 
rate versus /ihaio and function Mgai(Mhaio)- These shapes 
evolve weakly with redshift, as such our results can be reason- 
ably understood with a simple instantaneous calculation. Tak- 
ing the z — known Mgai(Mhaio) alone, we can calculate the 
relative contribution to the differential growth rate of bulges 
at z = as a function of mass (essentially this amounts to 
populating a z = simulation with the observed HOD, and 
evolving it forward for some arbitrarily small amount of time, 
then calculating the relative importance of mergers as a func- 
tion of their mass ratio). Figure |9] compares the results from 
this simple procedure with an integration over merger history 
(of course, without the full merger history, we can only de- 
fine this in terms of the contribution to bulge assembly, not 
formation). There is little difference between the two. Be- 
cause at any given stellar mass the most important mergers 
are those that happened while the system was relatively near 
that mass (not mergers that happened when the system was 



much lower mass, since those by definition will contribute lit- 
tle to the present total mass of the system), the "memory" of 
early formation or growth at low masses is effectively erased. 
This makes our results robust to details of the model at low 
masses and/or high redshifts, where empirical constraints are 
more uncertain. 

5.2. Dependence on Galaxy Mass 

We have shown how the relative importance of different 
mass ratio mergers depends on mass in Figures [8] & |9l Fig- 
ure [TO] summarizes these results. We plot the median (/igai) 
(specifically the mass-weighted median /igai, corresponding to 
the merger mass ratio above which 50% of the mass in bulges 
was assembled; i.e. where /buige(> /^gai) 0.5 in Figure |9|l 
as a function of stellar mass, at z = 0. We also plot the cor- 
responding ±1(7 and 10 — 90% ranges. As demonstrated in 
Figure|9l similar results are obtained with an "instantaneous" 
calculation; in §|6l we show similar results varying a number 
of choices in the model. 

As seen before, major mergers dominate near ^ 10'" — 3 x 
10" M0, with minor mergers increasingly important at lower 
and higher masses. (In terms of the initial bulge formation, 
rather than assembly, the prediction would be the same but 
without the "turnover" at high masses - i.e. asymptoting to a 
constant (/igai) ^ 0.5 at high masses.) Most of the variance 
comes from differences in merger histories at fixed mass. At 
all masses, the range of contributing mergers is quite large 
- there is always a non-negligible contribution from minor 
mergers with mass ratios ^0.1—0.3. 

We stress that the turnover at high masses does not come be- 
cause of fewer major mergers. In fact, we have shown explic- 
itly that the number of major mergers in the primary history 
increases monotonically with galaxy mass. Rather, at high 
masses, the number of minor mergers increases even faster. 
Thus the relative importance of minor mergers is enhanced 
at the highest masses (picture for example the growth of a 
BCG via accretion of many satellites in a cluster). This is why 
the turnover would not appear if we made Figure [TO] in terms 
of the mass ratios important for bulge formation, instead of 
bulge assembly. So we expect these systems to be more bulge 
dominated as compared to ^ galaxies, but with interesting 
second-order effects in e.g. their kinematics and light profile 
shapes that indicate the role of many minor mergers in their 
recent history. 

5.3. Dependence on Bulge-to-Disk Ratios 

Figure [TT] examines how the distribution of contributing 
/igai (in terms of bulge assembly) scales with the bulge-to- 
total stellar mass ratio B/T of galaxies. At fixed mass, galax- 
ies with higher B/T are formed in preferentially more major 
mergers, and the trend is similar at all masses. This is the nat- 
ural expectation: a more major merger yields a system with 
higher B/T. Because /igai dA^merger/d log /igai is not quite flat 
in /igai (rising to larger fS), and ~ 1 significant mergers are 
expected since z ~ 2 (i.e. at times when the galaxy is near is 
present mass), the local B/T will be dominated by the largest 
merger the system has experienced in recent times. Many ob- 
jects all have some amount of bulge built by /igai 0.1 merg- 
ers - the question is which will have larger mergers that con- 
vert more mass to bulge. 

Figure [12] summarizes these results, showing the (mass- 
weighted) median merger mass ratio (/igai) contributing to 
bulge formation as a function of B/T at different stellar 
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Fig. 10. — The relative importance of different mass ratios to bulge assembly as a function of z = stellar mass. We plot the (mass-weighted) median ^j^i 
contributing to the assembly of bulge mass density (defined as fi„.^[ where /buisc[> Mgall = 0-5) as a function of stellar mass, along with the interquartile range 
and 10 — 90% range (lines bracketing different ranges as labeled; horizontal dashed line denotes the standard "major merger" definition fi > 1/3). Error bars 
show the range resulting from variations to the model (see §|6). Recall, the number of major mergers increases monotonically with mass; the decrease here at 
high mass is simply because the number of minor mergers increases yet more rapidly. 



masses. At all masses, even masses where the global bulge 
population is predominantly formed in minor mergers, galax- 
ies that are bulge-dominated (the E/SO population) are pre- 
dominantly assembled (and formed) in major mergers. In 
principle ten 1:10 mergers in a short time will form as much 
bulge as a single 1:1 merger. However, 1:10 mergers are 
not ten times more common, and as such are not an impor- 
tant or efficient channel for the formation of bulge-dominated 
galaxies. Recall that the average galaxy still experiences 
only ~ 1 min or 1:1 merg ers since z ^ 2 (see Figure |3]& 
IStewart et all (l2008l l2009a l)): the case of ten 1:10 mergers is 
then a ~ 5 — lOcr outlier. Moreover, even if a system has sev- 
eral such mergers, they will be spaced widely in time (they 
essentially never occur simultaneously), so the galaxy disk 
will re-grow, reducing B/T after each and offsetting the bulge 
growth from mergers. In contrast, ~ 1/2 of all galaxies un- 
dergo a single 1:3 merger since z^2; these will immediately 
form a large B/T system. In short, minor mergers are not so 
much more common than major mergers as to dominate the 
formation of high B/T systems. 

This is also important for reproducing the existence of disks 
(especially "bulge-less" disks) - if minor mergers were so 
common as to dominate the formation of high B/T systems (if 
e.g. half of ^ L^: galaxies had formed through the channel of 
^ 10 rapid 1:10 mergers), it would be correspondingly much 
more rare for a system to have undergone very/ew 1:10 merg- 
ers, necessary to explain the existence of at least some signifi- 
cant number of low-mass systems with B/T <Q.\. Moreover, 
in practice any system with such an extreme merger history is 
likely to have also experienced an enhanced major merger rate 



- so the major mergers will still dominate bulge formation (it 
is unlikely to contrive an environment with so many ~ 1 : 10 
mergers in a short time and no major mergers). 

On the other hand, this implies that minor mergers do dom- 
inate the formation of bulges in low B/T galaxies (Sb/Sc/Sd 
galaxies). This is for the same reason - most galaxies have 
experienced ^1 1:10 merger in recent times (z < 2), whereas 
only some fraction have undergone more major mergers. The 
"traditional" scenario for bulge formation - early formation in 
a major merger, followed by subsequent re-growth of a disk 
by new cooling - is only responsible for a small fraction of 
the mass density in disk-dominated B/T < 0.2 systems. It is 
very rare that a system would have such an early major merger 
but then not have a later ^1:10 merger in the Hubble time 
required to grow the galaxy by a factor ~ 10 in mass. 

These results are independent of all our model variations 
(§ ISJ, so long as we ensure that we reproduce a reason- 
able match to the observed HOD and halo merger rates. In 
fact, these conclusions appear to be quite general, similar to 
those found from other models that ad opt different models for 
bulge form ation in mergers (see e.g. iKhochfar & Silkll2008l: 
IWeinzirl et al. 2009). 

Figures[TT]&[T2]also show how, at fixed B/T, the contribut- 
ing mass ratio distribution depends on stellar mass. At fixed 
B/T, the trend with mass is much weaker than seen compar- 
ing all galaxies as a function of mass (Figure [TOll. Moreover, 
the trend at fixed B/T appears to have an opposite sense: low- 
mass galaxies require higher mass ratio /igai mergers to reach 
the same B/T . This is priarily a consequence of the depen- 
dence of gas fraction on stellar mass and the effects of gas 
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Fig. 1 1 . — Left: Contribution of mergers of different fj,g^\ to bulge assembly 
(as Figure|9] differential version), as afunction of bulge-to-total ratio B/T at 
fixed galaxy mass. Line type denotes the final B/T value (as labeled), and 
each panel shows galaxies of a different mass. At all masses, more bulge- 
dominated systems are formed by more major mergers. Ellipticals and SO's 
are dominated by major merger remnants; late-type disk bulges are preferen- 
tially formed in situ in minor mergers. Right: Same, as a function of galaxy 
stellar mass at fixed B/T (lines denote different galaxy masses, panels show 
results for systems with different final B/T). At fixed B/T, the residual de- 
pendence on mass is weak; low-mass galaxies are more gas-rich, so require 
more major mergers to reach the same B/T. 




Fig. 12. — Summary of Figure l 1 1 1 median mass ratio of mergers contiibut- 
ing to bulge growth (defined as in Figure [To) as a function of the final B/T 
and galaxy stellar mass (different lines). Roughly speaking, the correlation 
reflects the instantaneous scaling for a single merger, B/T oa (as ex- 
pected if galaxies grow in a manner such that the most recent mergers, at e.g. 
z < 2, are most important). 
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Fig . 13. — As Figure fm but showing the confiibution to bulge growth from 
mergers with different gas fractions /gas. Top: Integrated distribution (frac- 
tion of bulge mass formed in mergers with /gas above the given value; lines 
show results at different stellar masses, with the "all" line integrated over the 
entire bulge mass function). Bottom: Same, in differential form (contribution 
per unit /gas). We compare the median z = 2 gas fractions of disks of the 
same mass (diamonds, color corresponds to the mass as labeled). To low- 
est order, the gas fractions of bulge-forming mergers simply reflect the gas 
fractions of disks at the time of merger (z = 2 is just representative; the dis- 
tribution of merger times is broad). The bulge mass density is dominated by 
mergers with /gas ~ 0. 1 — 0.2, with a tail towards more dissipational mergers 
in lower-mass systems. 

on bulge formation (Figure |7]). A low-mass galaxy, being 
very gas-rich, might require a major merger to even get to 
B/T ^ 0.2 (if, say, /gas ^ 0.8) - so low-mass systems will 
require more major mergers. On the other hand, mergers in 
a massive ^ IO^^Mq system, being gas-poor (/gas < 0.05), 
will yield B/T > 0.2 for any mergers with /igai > 0.2; the 
^ IO'^Mq galaxies with low B/T (what few there are) must 
be those that had only minor mergers in the last few Gyr. 

Overall, however, we wish to stress that systems with large 
B/T at low masses and low B/T at high masses are rare - 
most low mass systems, having low B/T, have had relatively 
more contribution to their bulge growth from minor mergers, 
and most higher mass systems, having high B/T, have had 
increasing contributions from major mergers. 

5.4. Dependence on Galaxy Gas Fractions 

Figure [13] compares the contribution to bulge formation 
from mergers not as a function of mass ratio, but as a func- 
tion of the gas-richness of the merger, where /gas is here 
defined as the sum gas-richness just before the merger (= 
(Mgas,i +Mgas,2)/(A^gaii +A^gai,2))- In an integrated sense. 
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the most important gas fractions for bulge formation are /gas ^ 

0.1 —0.20 This agrees well with estimates from numeri- 
cal simulations of the gas fractions required to form realistic 
^ ellipticals (in terms of their profile shapes, effective radii 
and fundamental plane correlations, rotation and higher-order 
kinema tics, isop hotal shapes, triaxiality, and other properties; 
see e g.lNaab et al. 2006; Cox et al. 2006a; Je sseit et al.ll2007l 
l2009t iHopkiiiTet al., i2008ci.a. ,2009c Fl In a cosmological 
sense, it simply reflects the gas fractions of ^ disks, the 
progenitors of ^ ellipticals. We stress that this does not 
mean the bulges are made purely from this (relatively small) 
gas mass - rather, this represents the typical mass fraction 
formed in a central starburst in the bulge-forming merger; the 
majority of the bulge mass is formed via violent relaxation of 
the pre-merger disk stars. 

As a function of mass, the typical merger contributing to 
bulge formation is more gas-rich at low masses. But as shown 
in Figure [13] this largely reflects the trend of gas fractions in 
late type or star-forming galaxies as a function of mass. At 
a given mass, in particular at the lowest masses where gas 
fractions can be sufficiently high as to significantly suppress 
bulge formation, there is a weak tendency for the dominant 
mergers contributing to bulge formation to be less gas-rich 
(since such mergers will, for the same mass ratio, form more 
bulge). However, the effect is not large. 

An important check of this is that it reproduce the "dissipa- 
tional" mass fraction in observed ellipticals, as a function of 
mass. This is the mass fraction of the spheroid formed in a 
dissipational starburst, rather than violently relaxed from the 
progenitor stellar disks. Being compact, this component is 
the primary element that determines the effec tive radii, pro- 
file shape, and ellipticity of a merger remnant. 'Hopki ns et alJ 
(|2009a e, 2008a) develop and test an empirical method to es- 
timate the dissipational mass fraction in observed local el- 
lipticals, and apply this to a wide range of observed ellip- 
ticals waA_j_combinationof HS T and ground-base d data 
from iKormendv et all (l2009h and iLauer etalJ ( l2007h . As 
shown therein, resolved stellar population studies yield sup- 
porting conclusions. Figure [T4l shows these empirically in- 
ferred dissipational fractions as a function of mass, and com- 
pares the predictions from the models here. The agreement 
is reasonable. Similar conclusions are reached even by mod- 
els with significantly d ifferent bulge formation prescriptions 
jKhochfar & SiUd '2006'). Note that the observed systems here 
are all classical bulges, appropriate for comparison to our pre- 
dictions. 

We also compare observed disk gas fractions. To lowest 
order, the dissipational fractions simply trace these gas frac- 
tions, but at low masses, the predicted and observed dissipa- 
tional fractions asymptote to a maximum ^ 0.3 — 0.4. This is 
because angular momentum loss in the gas becomes less effi- 
cient at these high gas fractions; if the fraction of gas losing 
angular momentum scales as adopted here, then the dissipa- 
tional fraction of the bulge formed from disks with gas frac- 

Recall, this is the gas fraction at the time of the merger, and can be 
different from the "initial" gas fraction at the beginning of an interaction, 
depending on e.g. the efficiency of star formation and stellar feedback. For 
example, in hydrodynamic simulations of idealized mergers without ongoing 
continuous accretion, a gas fraction at the time of merger of ^ 0.1 —0.2 
corresponds to an "initial" gas fraction ~ 2 Gyr before merger of ^ 0.3 — 0.4. 

In fact, only mergers with these properties have been shown to yield a 
good match to these quantities: mergers with significantly less or more gas, 
as well as secular instabilities and dissipational collapse have been shown to 
yield remnants with properties unlike observed ellipticals. 
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Fig. 14. — Top: Dissipational fraction (mass fraction of bulges formed 
in starbursts from gas that has lost its angular momentum in mergers, rela- 
tive to the total - starburst plus violently relaxed former stellar disk - bulge 
mass) as a function of stellar mass (solid line is the predicted median; dotted 
the ±1 (T scatter). We compare to empirically inferred /dissipational from de- 
composition of high-resolution surface brightness profiles and kinematics of 
observed elUpticals, presented in Hopkins et al. 1 2009a e) with samples from 
iKormendv et aU 1,2009) and Lau er et aL (.200 7). Bottom: Predicted dissipa- 
tional fraction from above (solid line), and median gas fraction /gas of merg- 
ers contributing to bulge growth (from Figure [Til dashed line), compared to 
observed disk gas fractions (from FigurefT] points in the same style, with val- 
ues at z = and j = 2 in blue and red, respectively). Dissipational fractions 
reflect the gas fractions of progenitor disks, but with an asymptotic upper 
limit of /dissipational ^ 0-4 that reflects the suppression of angular momentum 
loss in very gas-rich mergers. 

tion /gas is not /gas, but ~ /gas/(l +/gas), 1-6. asymptoting 
to the values observed for all /gas ~ 0.5 — 0.9. 

5.5. Redshift Evolution: Can Mergers Account for the Mass 
Density in Bulges? 

At all redshifts, the distribution of /igai contributing to 
bulges is similar to that at z = 0, shown in Figure [15] The 
only significant evolution is that the "turnover" in (/tgai) at 
high masses for the mass ratios contributing to bulge assem- 
bly (Figure[TOll becomes less pronounced and moves to higher 
masses. Technically, this relates to how Mgai(Mhaio) (empiri- 
cally constrained) evolves, but physically it is simply under- 
stood: at higher redshifts, "dry" assembly is less important, 
so the assembly (/Xgai) increasingly resembles the formation 
(Mgai) (Figure ISJ. By z ^ 2, there is no difference - dry as- 
sembly is negligible, and all high-mass, high-B/T systems 
are preferentially formed in major mergers. 
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Fig. 15. — As Figure[lO] showing the (mass-weighted) median mass ratio 
^igai contributing to bulge assembly as a function of mass for samples at dif- 
ferent observed redshifts, z = — 2 (different hnes, as labeled; again the red 
horizontal line denotes the traditional major merger definition /i > 1 /3). The 
qualitative trends are similar; at high redshifts the break/turnover moves to 
higher masses. 




Redshift z 

Fig. 16. — Predicted integrated mass density in bulge-dominated galax- 
ies (B/r > 0.4; black solid line) as a function of redshift, compared to 
observations (points). Obser vations are from the mo rphologically-selected 
samples of Bell et al] <200l black x), IBundv et ali (12005, 2006, red cir- 
cles), Abraham et al. ( 2007, violet diamonds), and Daddi et al. (2005, cyan 
squ are), and color-selected samples of Labbe et al. (2005, green square), 
Ivan Dok kum et al. ( 2006, orange square), and iGrazian et alj ( i2007l magenta 
stars). W e compare the total stellar density obsei'ved (green dot-dashed line; 
IHopkins & Beacom 2006). Increasing gas fractions and cooling rates at high 
redshift suppress the bulge mass density relative to the total stellar mass den- 
sity. The predicted number of mergers is sufficient to account for the z = 
and high-redshift evolution in the global bulge mass budget, but with factor 
~ 2 uncertainties. 

Given these predicted merger rates and BjT distributions, 
we also obtain a prediction for the mass density in bulge- 
dominated galaxies as a function of redshift. Figure [16] com- 
pares the redshift evolution of the bulge mass density to that 
observed. E3 The agreement is good: not only are there a 

Specifically, we plot the mass density in bulge-dominated galaxies, 
which is not the same as the absolute mass density in all bulges, but is closer 
to the observed quantity. At high redshifts z > 1.5 observed morphologies 
are ambiguous; we show the mass density in passively evolving red galaxies 
as a proxy. This may not be appropriate, but at z < 1 the two correspond 
well, and the compactness, size, and kinematics of the "passive" objects do 
appear d istinct from star-forming objects ( Kriek et al. 2006 ; Toft et al. 2003; 
ITriiiillo et al. 2007; Franx et al. 2008; Genzel et al. 2008). Also, the observa- 
tions in some cases do not distinguish "classical" and "pseudo" bulges; but as 
we discuss in §|7] the latter where measured appear to contribute only ~ 10% 



sufficient number of major mergers to account for the ob- 
served merger fractions, but also to account for the observed 
buildup of the bulge population with redshift. This should 
not be surprising, given the agr eement with observed merger 
fractio ns demonstrated in § 14.21 above: Hopkins et al.' ("200731 
|2008bj) . Hundyetal. (2009), and Bell et al. (2006a) have 
demonstrated that observed major merger fractions are suffi- 
cient, within a factor ^ 2, to account for the observed growth 
of the bulge population over the same redshift interval (given 
the observable lifetime calibrations that we adopt in § 14.21 
note that some of these works use different merger timescale 
estimates and reach different conclusions, but they are consis- 
tent using a uniform, simulation-calibrated timescale). Simi- 
lar results as a function of galaxy morpholog y are suggested 
by local observations (e.g. iDarg et al.ll2"009al) . Likewise the 
agreement between the predicted integrated number of merg- 
ers an d various observation al estimates suggests th is is plau- 
sible dde Ravel et al.1 l2009t IConselice et al. 20091: iLin et alj 
l2008h . 

For a more detailed com parison, as a function of e.g. galaxy 
stellar mass, we refer to IHopkins et al.l (l200 9g). who use 
the same merger rates as modeled here to predict e.g. the 
morphology (B/r)-mass relation and bulge/disk mass func- 
tions as a function of redshift. Provided proper account of 
galaxy gas fr actions is taken, good agreement is obtained. 
IStewart et al.l (l2009b) perform a similar calculation (with a 
basic criteria for bulge formation), with merger rates in close 
agreement as a function of galaxy mass to those measured ob- 
servationally in Bundv et al. (2009), and obtain similar good 
agreement with the bulge mass function as a function of red- 
shift. They actually find that bulge mass is somewhat over- 
prod uced, without a ccounti ng f or the role of g as-rich merg- 
ers. IHopkins et al] (l2008dl) and IHopkins et all (|2008b) con- 
sider a range of model parameter space (with several of the 
specific model variations discussed in § |6]l and perform sim- 
ilar calculations; they explicitly show the predicted spheroid 
mass function and mass fraction as a function of stellar mass, 
halo mass, environment, and redshift, for the different mod- 
els considered. They likewise conclude that, for all the model 
variations considered (with scatter in merger rates as a func- 
tion of mass similar to that discussed below), good agreement 
with the mass function and mass density of classical bulges, 
at masses > lO'^M©, is obtained. At lower masses, however, 
uncertainties grow rapidly. 

Relative to the total stellar mass density, the mass density 
in bulge-dominated galax ies decrea ses with redshift. This is 
discussed in detail in Hopkins et al.' (2009g), but the reason is 
simply that at higher redshifts, higher gas fractions suppress 
bulges (and the suppression moves to higher mass, relative 
to the galaxy mass function break. This trend agrees with 
that observed and is not trivial (models neglecting the impor- 
tance of gas-richness in affecting bulge formation efficiency 
in mergers, for example, may predict the opposite). 

5.6. Analytic Fits 

It is convenient to fit the distribution of merger mass ratios 
contributing to bulge formation at a given mass. The average 
^gai contributing to bulge assembly, i.e. (/-igai) where /bulge (> 
Mgai) = 0.5, as a function of z = galaxy stellar mass 
(Figure [Toll can be approximated as 



^"^^"^ (M,/1O11M0)-"-5 + (M,/1O11Mq)"-S ■ 

of the mass density. 
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If instead the (/igai) contribution to bulge formation is desired, 
a similar formula applies, but with a weaker turnover at high 
M*, i.e. 



1 



gal/ 



(M*/1011Mq)-"-5 + (M,/1011Mq)"-2 ■ ^^^^ 

In greater detail. Figures [TTIfTSl demonstrate that the typical 
(Meal) contributing to bulge assembly depends on the bulge 
mass fraction B/T at a given mass. As a bivariate function of 
B/T and mass, (/igai) can be approximated by: 



1 + (M,/1011Mq)0- 



(13) 



Finally, knowing (/^gai), we find that the complete distribu- 
tion /buige(> Mgai) (e-g- Figure|8]l can be simply approximated 
by 



/bulge (> A^gal) - (1 -Mgal)'' 



(14) 



where (since /buige(> (Mgai)) = 0.5), the slope 7 is trivially 
related to (/igai) as: 



-hi2 



ln(l - (^gal)) 



(15) 



6. HOW ROBUST ARE THESE CONCLUSIONS? A COMPARISON OF 

MODELS 

The relative importance of e.g. minor and major mergers in 
bulge assembly owes to the combination of reasonably well- 
determined halo merger rates and halo occupation statistics. 
Nevertheless, there are still uncertainties in this approach. 
We therefore examine how robust the conclusions here are 
to a variety of possible model differences. A much more 
detailed investigation of e.g. differences in predicted merger 
rates between semi-empirical models, semi-analytic models, 
and simulations will be the subject of a companion paper 
dHopkins et al.l l2009f>. Here, we wish to examine differences 
arising within the semi-empirical framework. 

Figures fTTlfTSl compare our "default" model with a number 
of alternatives. For each alternative, we show the the merger 
rate as a function of redshift, and the integrated contribution 
of different mass ratios /bulge (> /^gai) to bulge assembly (inte- 
grated over all bulge masses). For the merger rates, we also 
compare with the observational constraints: the dotted region 
in the Figure shows the approximate ±1 ct allowed range from 
the observational compilation in Figures |4l|5] (fitting a piece- 
wise broken power-law to the observations), given the merger 
lifetime calibrations discussed in § 14.21 

The model variations we consider include: 

6.L Halo Occupation Models 

Here, we consider an otherwise identical model, but adopt 
a different set of halo occupation constraints to determine 
M^, (Mhaio) (for now, we keep Mgas (M» ) fixed, but varying that 
is very similar to va ryingM,,(Mhain))- Fifst, our default model, 
using the fits from IConrov & Wechsled (120091) . Second, the 
first to M, (Mhaio) and its scatter for central and satellite galax- 
ies fr om the observed SDSS clustering at z = (Wan g et aLl 
12006*): here, we simply adopt the z = fit at all redshifts - we 
do not allow for evolution. Third, assigning galaxies to ha- 
los and subhalos based on a monotonic rank-ordering method 
(see lConrov et al.l2006l) . fitting to the redshift- dependent stel- 
lar mass function from iFontana et al.l (12006b. Fourth, the 
same, w ith the mass functions from iPerez-Gonzalez et al.l 
(I2008bl) . 



We have also consider ed various fits dir e ctly taken from 
other sources, inclu ding lYan et all (I200I: ICqo rav (20051 
2006"); 'Conrov et al.' (2006'. '2007h: IZheng et al.r(l2007) and 
Perez-Gonzalez et al. ( 2008a); these he within the range 
shown in Figure [17] Using the HOD predicted by semi- 
analytic models, at least for central galaxies, also appears 
to give similar results (we hav e com p ared the results in 
Croton et all 120061: iBower et al.1 120061: Ide Lucia & Blaizoti 



2007h : given that these models are constrained to match the 



observed stellar mass function, this appears to be sufficient 
for the level of convergence shown. 

6.2. Merger Timescales 

In our "default" model, we assume a delay between halo- 
halo and galaxy-galaxy mergers, given b y the dynamical fric- 
tion ti me calibrated to simulations in iBoylan-Kolchin et al.l 
(l2008h . We now allow this to var y according to five dif - 
ferent scalings, described in detail in lHopkins et al.l (l2008dl) . 

(a) Dynamical Friction: the traditional dynamical friction 
time, using the calib r ation from numer ical simulat ions in 
IBoylan-Kolchin et al.1 (l2008h (see also iJiang et al.1 I2OO8,) . 

(b) Group Capture: the characteristic timescale for pair- 
pair gravitational capture in group enviro nments, ca hbrated 
to siniulations following iMamonI (12006b (see also IWTiitd 
119761: iMakino & H"u3 11997') (c) Angular Momentum Cap- 
ture: as group capture, but considering capture in angu- 
lar momentum space rath er than gravitationally, following 
iBinney & Tremain3(ll987l) . (d) Gravitational Cross-Sections: 
similar to the group capture timescale, this is the timescale 
for gravitational capture between passages in e.g. loose 
group or field environments, calibrated from simulations in 
Krivi tskv & KontorovichI (11997 ). (e) No Delay: simply as- 
suming galaxy-galaxy mergers occur when their parent halo- 
halo mergers do. 

Although the dynamical friction time is most commonly 
adopted in e.g. semi-analytic models, each of these timescales 
depends on certain assumptions and is relevant in different 
regimes. A dynamical friction time is appropriate for a small, 
dense satellite at large radii; it becomes less so at small radii. 
A group capture or gravitational capture cross section, on 
the other hand, is appropriate for collisions in small groups 
or field environments where "inspiral" is not well-defined. 
The an gular momentum calibration from Binnev & Tremain^ 
(11987b is more appropriate for satellite-satellite mergers. In 
any case, we see the choice makes little difference to our con- 
clusions. And the range between these choices is generally 
much larger than other, more subtle details (e.g. allowing for 
continuous satellite mass loss in inspiral, or resonant baryonic 
effects that speed up the coalescence). The reason is simply 
that we are largely focused on fairly major mergers, for which 
the merger time is short relative to the Hubble time. The rate- 
limiting step is the accretion of such a companion, not the 
inspiral time. 

6.3. SubHalo Mass Functions/Substructure-Based 
Methodologies 

Instead of using a halo-halo merger rate with some "de- 
lay" applied, we can attempt to follow subhalos directly af- 
ter the halo-halo merger, and define the galaxy-galaxy merger 
when the subhalos are fully merged/destroyed. This will self- 
consistently allow for some distribution of merger times ow- 
ing to e.g. a range of orbital parameters, and will include 
satellite-satelhte mergers (neglected in our default model). 
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Fig . 17. — Top: Predicted major merger rate (/Jgui > 1/3), varying the assumptions in the models (different lines). Black line is the "default" model assumption. 
Red dotted range is the approximate range allowed by observations (the compiled points in Table [T]& Figure |4l for clarity, we show the range of these points 
rather than each individual measurement). Bottom: Corresponding cumulative contribution of different mass ratio mergers to the assembly of the bulge mass 
density (as Figure[8] integrated over all bulge masses). Left: Changing the halo occupation constraints: the three cases from Figure[3]are shown, together with 
adopting the ; = SDSS fits from Wang et al. ( 2006) at all redshifts. Center: Changing the merger timescale between halo-halo and galaxy-galaxy merger: using 
dynamical friction times (calibrated in Bovlan-Kolchin et al. 2008); angular momentum-space capture cross sections from Binnev & Tremaine 1 1987); coUisional 
group capture cross sections from Mamon (2006); gravitational capture cross-sections for field/small group crossings from Krivitskv & Kontoroviclj U997i) : or 
no delay. Right: Tracking subhalos to assign merger times or using subhalo mass functions instead of halo merger trees as a starting point: usin^ the subhalo 
merger trees from cosmological simulations in , Stewart et alj|2009ai): adopting the subh alo mass functions from simulations in lKravtsov et aU t2004l) : or the same 
from extended Press-Schechter trees constructed foUowing lvan den Bosch et alj 4200^ . 



which IWetzel et al.l ( l2009ah show can be important at the 
~ 10 — 20% level independent of halo mass. 

Here, we compare our default model to those obtained 
tracking the haloH-subhalo populations in cosmological sim- 
ulations from Stewart et al. (2009 a) (popu lating subhalos ac- 
cording to our default HOD). IWetzel et a l. (2009a b) also an- 
alyze subhalo merger rates, with a different methodology. 
They reach similar conclusions, but with a systematic factor 
^ 1.5 — 2 lower merger rate. As they discuss, this is quite 
sensitive to how one defines e.g. subhalo versus friends-of- 
friends group masses; some of those choices of definition will 
be "normalized out" by the appropriate HOD (renormalized 
for whatever subhalo populations are identified in a simula- 
tion so as to reproduce the observed clustering and mass func- 
tions), but it also reflects inherent physical uncertainties in the 
instantaneous mass and time of subhalo merger. 

We also compare with the results usi ng the differ- 
ent subhalo-based methodology described in Hopki ns et alj 
(l2008 d) (essentially, beginning from the subhalo mass func- 
tion constructed from cosmological simulations and evolv- 
ing this forward in short time intervals after populating it, 
at each time, according to the HOD constraints). We com- 
pare two different constructions of the subhalo mass func- 
tio ns: that from cosmol ogical dar k-matter only simula tions 
in iKravtsov et"!!] (|2004 (see also IZentner et al] l2005h and 



that from the extended Press-Schechter formalism coupled 
to basic p rescriptions for subh a lo dyn amical evolution, de- 
scribed in Ivan den Bosch et al] (|2005|). Alternat i ve subhalo 
mass functions from e.g. iDe Lucia et alj (120041) : iGao et all 
(2004.) ; .Nurmi e t al. (2006!) are consistent. 

6.4. Halo Merger Rates 

We can next vary the halo-halo merger rate s adopted. 
Our de fault model uses the merger rates in Fak houri & Mai 
( 12008b) , calibrated fr om the Millenium dark-matte r only cos- 
mological simulation dSpringel et al.ll20055l2006h . Since we 
use the full history, this is equivalent to the "per progenitor" 
merger rates defined from the same simulation in Genel et al] 
(2008). An alternative dark-matter simulation, of comparable 
resolution, with halo merger r ates determined using a different 
methodolog y, is described i n lStewart et alj (l2009ah. Another 
is found in Gottlober et al] (1200 Ih (see also lKravtsov et alj 
120041; IZentner et alj|2005h ; thev quantify the fit separately to 
field, group, and cluster environments. 

We can also compare with the merger rates from 
iMaller et al. ( 2006), determined from cosmological hydrody- 
namic simulations. Although it is well known that, without 
proper implementations of feedback from various sources, 
cosmological hydrodynamic simulations yield galaxies that 
suffer from overcooling (and do not reproduce the observed 
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Fig. 18. — Figure [TtI continued. Left: Clianging the lialo-halo merger trees: our default choice from the Millenium simulation analyzed in 'Fakhouri & Ma| 
(2p08b); numerical trees from an alternative high-resolution cosmological simulation (Stewart et al. 2008, 2009a); cosmological DM only simulations of field, 
group, and cluster environments (as labeled; Got tliiber et al.l l2001); or merger rates from cosmological SPH simulations (tracking galaxy-galaxy mergers but 
still re-populating them appropriately for the observed HOD; Mailer et al. 2006, upper and lower correspond to their high and medium-mass primary sample, 
respectively). Center: Changing the cosmology; our default WMAP5 (n,„,as)=(0.n, 0.81) cosmology, versus a WMAPl (0.27, 0.84), WMAP3 (0.27, 0.77), 
and "concordance" (0.3, 0.9) cosmology. 



halo occupation statistics), the galaxies in these simulations 
can still serve as "tracers" of halos and subhalos. This pro- 
vides a means to avoid the considerable ambiguities in defin- 
ing a halo merger (moreover in considering the delay between 
halo-halo and galaxy-galaxy mergers). Although the galaxy 
masses may not be correct, they are still tracers of where in 
the halo real galaxies should be, and therefore can be used 
to measure the halo merger rate. We do so by recalculating 
their merger rates after re-populating the galaxies appropri- 
ately (essentially renormalizing their predicted mass function 
to match that observed). 

6.5. Cosmological Parameters 

We can also vary the cosmological parameters 
and see if this makes a significant difference to 
our conclusions. We consider four sets of cosmo- 
logical parameters: a "concordance" model with 
(r2M,f^A,/?,cr8,n5)=(0-3,0.7,0.7,0.9, 1.0), the WMAPl 
(0.27, 0.73, 0.71, 0.84, 0.96) results of Spergel et al. (2003), 
WMAP3 (0.268,0.732,0.704,0.776,0.947) (Soergel et al.i 
l2007h. and WMA P5 (0.274, 0.726, 0.705, 0.812, 0.96) 
(iKomatsu e t al.ll2009h . It is prohibitively expensive to re-run 
the simulations for each case, and moreover the qualitative 
behavior is not expected to change (seen in e.g. lower- 
resolution dark-matter simulations). We simply renormalize 
the halo masses at all times to match the halo mass function 
and accretion history appropriate for the revised cosmological 
parameters (see e.g. .Neistein et al., i2006i) - the dominant 



effect is the predicted halo mass function shifting to higher 
masses with larger ag. However, because we use a halo 
occupation-based approach, the model is re-normalized to 
yield the same observed galaxy mass function and clus tering, 
so these differences are largely normalized out. El ahi et al.l 
(2008) show that the quantity of greatest importance for our 
conclusions, the normalized substructure mass function or 
(equivalently) dimensionless merger rate (mergers per halo 
per Hubble time per unit mass ratio) is almost completely 
independent of cosmological parameters including e.g. the 
power spectrum shape and amplitude over the range of 
variations here (not until one goes to much larger effective 
ris ~ 3 does one see this function change shape). 

6.6. Bulge Formation Prescriptions 

We can also consider variations in the physical prescription 
by which bulge mass is formed in mergers. Obviously this 
will not change the merger rates, but it could change the rela- 
tive importance of mergers of different mass ratios. However, 
we are tightly constrained by the results of A^-body simula- 
tions; since the physics determining gas angular momentum 
loss and violent relaxation are predominantly gravitational, 
there is little uncertainty in how much bulge should be formed 
in a given merger (given the appropriately normalized initial 
conditions of interest). Still, there are some differences in fit- 
ted prescriptions: we have re-calculated the results from our 
default model according to the approximate results of simu- 
lations from .Naab & Burkert (.2003.) and .Naab et al.. (.2007,) . 
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as well as lBournaud et al.1 (120051) and ldi Matteo et all (I2007I) . 
W e have also us e d the fi ts to the same s uite of simulations 
in 'Hopkins et al.' ('2009b') as presented in ICox et all (l2006bb 
and Cox et al. ( 2008). Note that the results in several of these 
works do not necessarily include a complete survey of param- 
eters such as e.g. mass ratio, orbital parameters, and gas frac- 
tion; where not given we interpolate between the results pre - 
sented based on the model outlined in' Hopkins et al.l (l2009bl) . 
In any case, the differences in quantities such as the abso- 
lute bulge mass (especially in gas-rich, low mass systems) and 
dissipational fractions (fraction of mass formed in starbursts, 
rather than violently relaxed from stellar disks) of ellipticals 
can be non-negligible, but the relative contribution of major 
and minor mergers is almost identical. 

This will be true, it turns out, in any model where the 
amount of bulge formed in a given merger scales roughly in 
linear fashion with the mass ratio. As such, even highly sim- 
plified models which ignore the role of gas fraction and orbital 
parameters, and/or only violently relax the primary in major 
mergers (but do destroy the secondary in minor mergers), will 
still obtain the same qualitative features in /bulge (> Mgai); see 
e.g. lKhochfar & SiUd (120061) . 

6.7. Combinations of the Above: Typical "Scatter" 

We have considered various permutations of the above 
models, amounting to ~ 700 total models; our conclusions 
are robust to these combinations. The interquartile range be- 
tween this sampling of models lies within the observationally 
allowed range in terms of the merger rate, and yields very lit- 
de scatter in /bulge (> A*gai)- 

To the extent that comparison of these models can be con- 
sidered "scatter" or reflective of uncertainties in the theoreti- 
cal predictions, the corresponding typical "uncertainties" are 
as follows: around ~ and at slightly higher masses, uncer- 
tainties are small - a factor of ^ 1 .5 in merger rates (at z < 2; 
uncertainties grow at higher redshifts as in Figure |3]l and 
smaller in /buige(> Mgai), with /buige(Aigai > 1/3) 60 - 80%. 
At factors less than a few higher and lower masses, these un- 
certainties increase to a factor ^ 2 in merger rate, and factor 
~ 1.5 in the importance of major versus minor mergers. At 
much lower masses (^ lO^M©), uncertainties in both grow 
rapidly - here, the halo occupation statistics are not strongly 
constrained. Moreover, /Xgai is the baryonic (not just stellar) 
mass ratio, and systems at these masses are increasingly gas- 
dominated, so uncertainties in Mgas(MH.) can strongly affect 
our predictions. 

7. DISCUSSION 
7.1. Approach 

We have used an extensive set of models to examine galaxy- 
galaxy mergers and to identify robust predictions for the rel- 
ative importance of mergers of different mass ratios for bulge 
formation. Although halo-halo merger rates have been rela- 
tively well-understood, mapping halo-halo mergers to galaxy- 
galaxy mergers is not trivial. There can be significantly 
more or fewer major or minor galaxy-galaxy mergers, rela- 
tive to halo-halo mergers; likewise, bulge growth can be dom- 
inated by preferentially more major or minor mergers than the 
growth of the host halo. 

However, there is hope. Numerical simulations are con- 
verging in predicting how the efficiency of bulge formation 
scales with merger mass ratio (and what the "appropriate" 
mass ratio to use in these calculations should be), giving a 



straightforward set of predictions for how much bulge should 
be formed in a given galaxy-galaxy encounter To lowest or- 
der, the amount of bulge formed scales linearly in the merger 
mass ratio, close to the maximal efficiency possible for minor 
mergers (Hopkins et al. 2009b). 

Meanwhile, observations are converging on relatively tight 
constraints on halo occupation models: namely, the stellar and 
gas mass of the average galaxy hosted by a halo/subhalo of a 
given mass. The correlation between galaxy stellar mass and 
halo mass is monotonic and, to lowest order, amounts to a 
simple matched ra nk-or dering of the two, with small scatter 
(e. g. rConroy^etanBooa) . 

7.2. Conclusions 

This convergence makes the time ripe to examine the 
consequences of galaxy-galaxy mergers on bulge forma- 
tion. To good approximation, the salient features of the 
merger rate distribution can be captured by convolving the 
theoretically determined halo-halo merger rate with the 
empirically determined halo occupation statistics. Given this 
simple, well-constrained approach, there are some robust 
predictions that are insensitive to most if not all model details: 

(1) Major-merger f/igai > 1/3) remnants dominate the 
integrated mass density of merger/interaction-induced 
bulges at all redshifts (Figures [SirTTTl. Minor mergers 
(1/10 < /igai < 1/3) do contribute a significant, albeit not 
dominant, fraction 30%) to the assembly of the total mass 
density. More minor mergers /itgai < 1/10 are not important 
(contributing < 5 - 10%). 

(2) This statement is significantly mass-dependent (Fig- 
ures |9][T0ll. Although the relative major/minor contribution to 
halos is nearly mass-independent, the mass-dependent HOD 
shape leads to a galaxy mass dependence: major mergers 
strongly dominate bulge production around ^ L, (where most 
of the bulge and stellar mass density of the Universe lies). At 
masses ^ L*, merger rates at all mass ratios are suppressed, 
and minor mergers are relatively more important (since 
Mgai cxM^jijjj at these masses, approximately, a 1:3 halo-halo 
merger becomes a 1 :9 galaxy-galaxy merger; so all mergers 
are shifted to lower mass ratio, suppressing the number at 
some fixed /igai and relatively suppressing major mergers). 
At higher masses > L,, merger rates are higher, and major 

1 12 

mergers are relatively more important (here Mgai cx M^^^^, 
so a 1:9 halo-halo merger becomes a 1:3 galaxy-galaxy 
merger; making mergers at all significant mass ratios more 
abundant and relatively increasing the importance of major 
mergers). This is true until >a few L», where minor mergers 
again become relatively more important owing to the rapid 
dropoff in the number of "major" companions (equivalently, 
since most of the galaxy mass is concentrated near ^ L,, 
most of the incoming mass density is weighted near this 
region as well, which is a major merger at < a few L* and 
minor merger above). These trends are quite general, and 
the relative "increased weight" of major mergers will occur 
wherever the "mass-to-light ratio" or formation efficiency 
A^haio/A^gai has a minimum {'^L^). 

(3) For massive galaxies, there is a difference between 
the mass ratios important for bulge formation (the mergers 
which initially converted some disk mass into bulge mass) 
and those important for bulge assembly (the mergers that 
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brought together the present-day bulge from any combination 
of pre-existing bulges and/or disks). At low masses, the two 
are equivalent (they are only different where "dry mergers" 
are significant). At high masses, the description above 
applies to assembly. Most mergers onto ^ systems are of 
already bulge-dominated galaxies (i.e. dry), systems which 
first turned their disk mass into bulge ("formed" the bulges) 
at progenitor masses near ~ L*, where major mergers are 
most efficient. As a consequence, most bulge mass at all 
> masses is formed in major mergers (albeit again with 
non-negligible contributions from minor mergers); however, 
bulges are assembled in increasingly minor (dry) mergers at 
larger masses. 

(4) The relative importance of major and minor mergers 
is also significantly morphology-dependent (Figures fTTI - 
fT2l i. Bulge-dominated (E/SO or B/T > 0.4) galaxies are 
preferentially formed in major mergers; later-type (Sb/c/d 
or B/T < 0.2) galaxies are preferentially formed in minor 
mergers. Despite the fact that simulations show that e.g. ten 
1:10 mergers can yield just as much bulge mass as one 1:1 
merger, cosmological models show that they are not ten times 
more common. Moreover, this many minor mergers would 
have to happen in a time much less than a Hubble time in 
order to successfully build a bulge-dominated galaxy, and 
this scenario is very unlikely (even at high redshift; minor 
merger rates may increase, but so do major merger rates). 
However, since just one or two 1:10 mergers are sufficient to 
account for a B/T < 0.2 bulge, this is a common formation 
channel for small bulges, in particular more common than a 
major merger at high redshifts that destroys the entire disk 
followed by a factor of ~ 10 subsequent disk re-growth (even 
if this occurred, it would take ^ a Hubble time, in which time 
a 1:10 merger would be very likely, and that merger would 
then dominate the final bulge mass). To lowest order, bulges 
of systems of bulge-to-total ratio B/T are characteristically 
formed in mergers of mass ratio /igai ^ B/T (EquationfTSTl. 

(5) Gas-richness, with high gas fractions /gas > 0.5, can 
dramatically suppress the global efficiency of bulge forma- 
tion (from mergers at all mass ratios), and the important 
implications of this for establishing the morphology-mass 
relation and allowing for a si gnificant popu l ation of low 
B/T systems is discussed in iHopkins et al.l (12009 d) and 
IStewart et al.l (|2009b). However, it does not affect the merger 
rate, and because the effects are not mass ratio-dependent, 
it does not significantly affect the relative importance of 
major/minor mergers. Because low-mass galaxies are typ- 
ically more gas-rich, they require somewhat more violent 
merger histories to reach the same B/T as a comparable 
high-mass galaxy (Figures [TTIfT2] ). To lowest order the 
gas fractions of progenitor galaxies that contribute to the 
observed bulge population, and the fraction of bulge mass 
formed dissipationally (by gas losing angular momentum in 
mergers and forming stars in concentrated nuclear starbursts) 
simply reflect the cosmological average gas fractions of 
progenitor disks corresponding to the same stellar mass and 
assembly times (Figures [T3lfT4l) . We have included the effects 
of gas on merger dynamics because it is known to be very 
important; however, given the above, our key prediction in 
this paper specifically would be not be dramatically changed 
if we ignored these effects throughout. 



(6) The predicted major merger rate (mergers per galaxy 
per Gyr) agrees well with observed merger fractions from 
z ^ — 2 (Figures [3]|5]) when one accounts for the observ- 
able merger timescale determined by applying the same 
observational methods directly to high-res olution galaxy- 
galaxy merger simulations (see e.g. i Lotz et al. 2008a). The 
corresponding rate is ^ 0.5 major galaxy-galaxy mergers 
per central galaxy per unit redshift (in these units, nearly 
redshift-independent), around ^ L*, and is mass-dependent as 
per conclusion (2): half to two-thirds of the ^ population 
has had a major merger since z ^ 2, but the fraction is a factor 
^3 — 5 lower at an order-of-magnitude lower stellar mass, 
and becomes one (with many galaxies having a couple such 
mergers) at a factor of a few higher stellar mass (Figures |2]& 
|6]). The merger rate as a function of galaxy-galaxy baryonic 
mass ratio /igai, redshift z, and primary stellar mass can be 
reasonably well fit by the simple functions in Equations ISlfTOl 

(7) Integrating over all mergers, the predicted merger rates 
yield good agreement with the growth of the mass density 
in bulge-dominated galaxies, from redshifts z = — 1.5 and 
(to the extent that color and morphology are correlated) the 
passive/red sequence population from redshifts z = — 4 
(Figure [T6l l. The typical uncertainties in both theory and 
observations are at the factor ^ 2 level; this is an interesting 
range discussed below. 

7.3. Robustness 

We have examined how these conclusions depend on a va- 
riety of choices, including the empirical HOD constraints, the 
global cosmology, halo merger rates, substructure tracking, 
and merger timescales, and find that they are robust (§|6l Fig- 
ures [T7]|l8] l. 

Varying the halo occupation model within the range al- 
lowed by observations including weak lensing, clustering, 
group dynamics, and abundance matching methods all yield 
similar conclusions. A very different halo occupation model, 
for example simply assuming Mgai oc Mhaio, would yield very 
different conclusions, but observational constraints are suffi- 
ciently tight that within the range allowed, resulting variations 
are small. Varying the cosmological parameters primarily af- 
fects the absolute abundance of halos of a given mass, not e.g. 
the shape of the merger rate function, and since the observed 
galaxy mass function is fixed, this difference is simply folded 
into the halo occupation model and does not change our con- 
clusions. 

Halo-halo merger rates, likewise, are sufficiently converged 
between different simulations such that they yield no large 
differences. However, a halo-halo merger is not a galaxy- 
galaxy merger Typically, one attempts to better approximate 
the latter by adopting either some merger timescale, repre- 
senting a delay corresponding to subhalo orbital decay before 
the galaxy-galaxy merger, or by following subhalos directly in 
high-resolution cosmological simulations. One can also use 
galaxy-galaxy mergers identified in cosmological hydrody- 
namic simulations, after re-normalizing their masses to agree 
with empirical constraints. Considering variations in each of 
these choices, we find that they have little effect on the shape 
of the merger rate function, hence little effect on the relative 
importance of major/minor mergers. Further, some apparent 
differences in the resulting merger rate owe purely to defini- 
tions, and are implicitly normalized out in the HOD. Never- 
theless, these different approaches do yield important system- 
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atic differences in the absolute merger rate, at the factor ^ 2 
level. 

Independent models adopting the halo occupation method- 
ology described here also obtain results in good agree- 
ment (see e.g. Zheng et al. 2007; Brown etaLj 120081 
iPerez-Gonzalez et al. 2008a; Stewart et al. 2009a^ How- 
ever, models that attempt to predict galaxy formation and 
merger rates in an a priori manner, such as e.g. cosmologi- 
cal hydrodynamic simulations and semi-analytic models, have 
reached various mixed conclusions - some in agreement with 
those here, some not, with significantly larger variation in 
the predicted galaxy- galaxy merger rates than the fact or ~ 2 
above (compare e.g. 'Weinzir l et al.l 12009 1: iParrv et al 
Mailer et al. 2006; Naab et al.l 120071; iGovernato et al. 



aLll200: 



2009 



2007 



Quo & Whita i2008; Somervil le et al.ll2008i) . The origins and 
impUcations of these difference s is ex amined in detail in a 
companion paper (Hopki ns et aDl2009a . 

However, the important point for our conclusions is that 
these methods are fundamentally different; they are not 
strictly tied to observed halo occupation constraints as are the 
models here. As such, they can yield very different predic- 
tions. For example, it is well-known that cosmological sim- 
ulations without feedback yield efficient star formation at all 
masses, such that the predicted halo occupation has a form 
more like Mgai cx Mhaio , and so galaxy-galaxy mergers will, 
in such a model (without re-normalizing masses) trivially re- 
flect halo-halo mergers. Some semi-analytic models, mean- 
while, have well-known discrepancies between predicted and 
observed populations of satellite galaxies, which propagate 
to the predicted merger rates. It is increasingly clear that 
these semi-analytic models have considerable difficulty re- 
producing the observations of the merger history (generally 
with the sense that the semi-analytic merger rates/fractions 
are lower than those observed; see e.g. Jogee et al. | 2008t 
B ertone & Co nselice 2009; Lopez-Saniuan et al. 2009iJ). In 
Hopkins et al.1 120091), we show that this indeed primarily 



owes to well-known mis-matches between the predicted satel- 
lite galaxy properties and galaxy mass functions in these mod- 
els, and those observed. As such, the semi-empirical approach 
can perform much better in explaining the observations (of 
course, the model here cannot predict satellite properties as 
can a semi-analytic model, but that is not its purpose). By 
adopting the halo occupation constraints directly from ob- 
servations, the semi-empirical model simply bypasses a ma- 
jor, well-known theoretical uncertainty in attempts to predict 
merger rates directly from semi-analytic models or cosmo- 
logical simulations. This does not mean the answer is "built 
in" implicitly somehow in our models here - what it does 
mean, however, is that the apparent discrepancy between ob- 
servations and other predictions of the merger rate owes not 
to some fundamental problem of ACDM, but rather to well- 
known difficulties in properly modeling the accretion and star- 
formation histories of galaxies. 

Carefully accounting for these distinctions, the different re- 
sults in various models can be understood. And in fact, de- 
spite differences in some quantitative predictions, many o f the 
qualitative conclusio ns are the same; iParrv et all (l2009h and 
IWeinzirl et alJ (l2009h demonstrate that different SAMs reach 
similar conclusions regarding how merger rates and the rel- 
ative importance of major versus minor mergers scale as a 
function of e.g. galaxy stellar mass and redshift. 

7.4. Outlook and Future Work 



Convergence in predicted merger rates among different the- 
oretical approaches, at the factor ~ 2 level or better, is a re- 
markable achievement. Unfortunately, obtaining greater con- 
vergence in theoretical predictions will be difficult. Applying 
constraints from empirical halo occupation approaches to e.g. 
cosmological simulations and semi-analytic models is impor- 
tant. Tighter observational constraints on the halo occupation 
distribution, in particular at low masses and at high redshifts, 
will allow semi-empirical models such as those in this paper 
to greatly extend the dynamic range of predictions (as well as 
putting strong constraints on a priori models for galaxy for- 
mation at these masses and redshifts). 

However, we have shown that these differences only ac- 
count for a fraction of the scatter in theoretical predictions 
- subtle details of how e.g. halos are defined and followed 
become important at this level. Moreover resolution lim- 
its and the absence of baryons in simulations (which does, 
at the level of uncertainty here, have potentially important 
effects o n the longevity and m erger timescales of subhalos; 
see e.g. [Weinberg et al.l 120081) limit all theoretical models. 
Ideally, high-resolution cosmological hydrodynamic simula- 
tions could form the basis for halo occupation models: avoid- 
ing ambiguity in identifying a galaxy-galaxy versus halo-halo 
merger by simply tracking the galaxies (even if their abso- 
lute masses are incorrect, and they need to be "repopulated" 
in post-processing). Although some steps have been made in 
this direction, it remains prohibitively expensive to simulate 
large volumes at the desired high resolution with gas physics. 

It is also unclear whether a merger rate alone is meaningful 
at an accuracy much better than a factor ^ 2. At this level, the 
question of e.g. the " proper" mass ratio becomes important 
(see e.g. IStewartI [20091) . What matters, in detail, for galaxy 
dynamics and the effect of a given merger is a combination of 
several quantities in the merger "mass ratio" - including stars, 
gas, and the tightly bound portion of the halo that has been 
robust to stripping; as such, the halo structure and history, as 
well as effects such as adiabatic contraction, become impor- 
tant. Moreover, at this level, the orbital parameters, galaxy 
gas fractions, and progenitor structure (relative bulge-to-disk 
ratios and disk scale lengths) become non-trivial corrections 
to the estimate of the effects "per merger." Without models 
for all of these details, a merger rate constrained to arbitrarily 
high accuracy does not necessarily translate to a bulge forma- 
tion model with accuracy better than a similar factor ~ 2. 

In the meantime, however, there is considerable room for 
improvement in the comparison of model predictions and ob- 
servations of the merger rate. The formal statistical errors 
in observed merger and close pair fractions are rapidly de- 
creasing; even including cosmic variance, such observations 
at z ^ — 1.2 are converging to better than a factor of ~ 2. 
However, as shown in Figure |4l simply putting all such ob- 
servations on equal footing yields order-of-magnitude scatter; 
similar uncertainties plague the conversion of these quantities 
to merger rates (and it is further unclear what the sensitiv- 
ity is to different mass ratio mergers). This is an area where 
considerable improvements can and should be made: most 
of the differences in observational estimates are attributable 
to different methodologies, observational depth, and selection 
effects. The conversion of some specific pair or morpholog- 
ically identified sample to a merger rate should be calibrated 
to suites of high-resolution A^-body simulations, specifically 
with mock observations matched to the exact selection and 
methodology adopted. 

Moreover, the merger rate is predicted to be a non-trivial 
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function of galaxy mass: many samples identifying merger 
fractions have ambiguous luminosity selection; what is ul- 
timately necessary are samples with well-defined stellar or 
baryonic mass selection. At the level of present data qual- 
ity and theoretical convergence, order-of-magnitude estimates 
of merger lifetimes and lack of such calibration represent the 
dominant uncertainty in comparisons. 

Impr ovements are being made in this area: iLotz et al] 
(l2008ah have calibrated the merger timescale for major 
pair samples of different separations and certain specific 
automated morphological selection criteria to mock obser- 
vations of_Jiigh;^esolutionhydrodynamic merger simula- 
tions. IConsehce et al.l (120091) adopt these and similar de- 
tailed calibrations to attempt to address consistency between 
merger populations identified with different methodologies. 
Upgee et al. (2008, 2009) attempt to calibrate their morpholog- 
ical selection criteria as a function of merger mass ratio. Var- 
ious works have attempted to quantify merger rates as a func- 
tion of stellar mas s, rather than in a pure magnitude- l imited 
sample (see e.g. 'Bell et al.' '200 6at IConselice et all 120081: 



IT992I; ISchweizeS fT992l; iBo urnaud et all llOOl iNaab et al 



2006; Naab & Trujillo 2006; Cox et all l2006at iJesseit e~ 



sample (see e.g . 'Bell et al.' 'ZUU bat IConselice et al.l UtJUbl: 
Bundv etal.1 12009; Lopez-Saniua n et al.1 l2009bl: iDarg et all 



2009b'). 

Together, these approaches will allow rigorous comparison 
of predicted and observed galaxy-galaxy merger rates as a 
function of galaxy stellar mass, redshift, and (ideally) mass 
ratio. Obviously, extension of observational constraints in 
any of this parameter space represents a valuable constraint 
on the models here. Using the calibrations above, we attempt 
such a comparison specific to different observational meth- 
ods (at least in terms of pair versus morphological fractions), 
and find good agreement between predicted and observed 
merger rates, and the integrated buildup of the bulge popu- 
lation. Considering the most well-constrained observations 
and well-calibrated conversions, we find agreement within a 
similar factor ~ 2 as that characteristic of the theoretical un- 
certainties. 

Far from implying that the problem is "solved," such a fac- 
tor of ^ 2 is of great interest. There is a large parameter 
space where predicted merger rates are consistent with ob- 
served merger/pair fractions as a function of mass and redshift 
and can be tuned to precisely account for the entire bulge mass 
budget of the Universe. However, allowing for the factor ^ 2 
uncertainty in one direction would lead to "too many" merg- 
ers, implying that mergers must be less efficient than cosmo- 
logically predicted: this might mean that real gas fractions are 
in fact higher than what we have modeled here, or that tidal 
destruction of satellites is efficient, even in the major merger 
regime, or that there is some problem in our understanding of 
halo occupation statistics or cosmological dark matter merger 
rates. 

On the other hand, allowing for the same factor of ^ 2 
variation in the opposite sense would imply that ^half the 
bulge mass density of the Universe could not be attributed to 
mergers as we understand them. This means that, within the 
present uncertainties, secular processes such as bar or disk in- 
stabilities might account for up to half of the bulge mass of 
the Universe. Since the uncertainties grow at low mass, the 
fraction could be even higher at lower masses. 

Independent observational tests can put complementary 
constraints on these possibilities. It must be emphasized, for 
example, that essentially all numerical studies of spheroid 
kinematics find that only mergers can reproduce the observed 
kinematic properties of observed elliptical galaxies and "clas- 
sical" bulges (.Hernquist .1989. .199Z .1993; .Barnes, [T9881 



l2007i) . These are, in general, the bulges whose formation 
history we predict here. Disk instabilities and secular evo- 
lution (e.g. bar instabilities, harassment, and other isolated 
modes) can indeed produce bulges, but these are "pseu- 
dobulges" ( Pfenniger 1984; Combes et al. 1990; Rahaetaj] 
T99?,'Kuii ken & Merrifield..l995; .O'Neill & Dubinski,2003T 
Athanassoulal l2005b . with clearlv distinct shapes, kinematics, 
structural properties, and colors from classical bulges (for a 
review, see Kormendy & Kennicutt 2004). 

Observations at present indicate that pseudobulges con- 
stitute only a small fra ction of the total mass density in 
spheroids (< 10% ; see lAllen et all I2OO6I; iBall et all I2OO6I; 
iDriver et"al] 12007 ^; they do, however, become a large frac- 
tion of the bulge population in small bulges in late-type 
ho sts (e.g. Sb/c, correspondin g to typic al M^ai S. 10"^ M0; 
see lCarollo et alii 998: Kormendy & Kennicutt2004. and ref- 
erences therein). However, this is not to say that secular 
processes cann ot, in principle, build some massive bulges 
(see e.g. Debat tista et al.l 12004 120061) . And it is not clear 
that mergers - specifically minor mergers with mass ratios 
Mgai ^ 1 /lO - cannot build pseudobulges, depending on e.g. 
the structural properties o f the secondary and orbital parame- 
ters of th e merger (see e.g. Gauthier et alJ2006l : lYounger et alJ 



l255^lEiIche-Moral et al.ll2008ir 

Improvements in theoretical constraints (from high- 
resolution simulations) on how bulges with different struc- 
tural properties are formed, combined with improved ob- 
servational constraints on the distribution of these structural 
properties, can constrain the role of secular processes at bet- 
ter than a factor ~ 2 level (at least at low redshifts) - a 
level at which theoretical models cannot yet uniquely pre- 
dict the importance of mergers. On the other hand, observa- 
tional constraints on the mass budget in extended galaxy ha- 
los, intra-group an d intra-cluster light can constrain satellite 
disruption (s ee e.g. 'Lin & Mohi^'2004* I Cypriano et~an 120061: 
iBrown et~aL 2008; Laganaetal. 2008), and observations of 
high redshift disk-nbulge systems that may represent recent 
re-forming or relaxing merger remnants can constrain the ef- 
ficiency of bulge formation in m ergers ( Hammer et al. 2005; 
|Zhen2 et al .1120051; iTruijUo & Pohlen 2005; Flores e t aL,2006l: 
Puech et al. 2007a b, 2008; Atkinson et al. 2007). Together, 
these improvements in observational constraints and theoreti- 
cal models have the potential to enable precision tests of mod- 
els for bulge formation in mergers, and allow a robust determi- 
nation of the relative roles of secular processes, minor merg- 
ers, and major mergers in galaxy formation, as a function of 
cosmic time and galaxy properties. 

In order to facilitate comparison with future observations, 
we have provided fitting functions to both the predicted 
merger rates as a function of galaxy mass and mass ratio, and 
to the relative contributions of these mergers to bulge forma- 
tion. However, for various applications, additional informa- 
tion is desired. We therefore make public a simple "merger 
rate calculator" cod43 which can be used to obtain the pre- 
dicted merger rates from the models as a function of e.g. 
galaxy mass, mass ratio, redshift, and galaxy gas fractions. 
The script can be used to determine merger rates as a func- 
tion of halo, stellar, or total galaxy baryonic masses, and can 
be used to restrict to e.g. gas-rich ("wet") or gas-poor ("dry") 



http : / / www . cf a . harvard . edu/~phopkins / Site/mergercalc . 
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mergers. It also allows for different choices with respect to 
e.g. the stellar mass functions used to normalize the HOD and 
Mgai(Mhaio) distribution used in the models here, and different 
cosmological parameters. As desired, it can output the merger 
rate per galaxy, the volumetric total merger rate (mergers per 
unit volume per unit time), or merger fractions with the appro- 
priate observable timescales used here as calibrated for pair or 
morphologically-selected samples. We note that, in the inter- 
est of running time and memory use, the script uses some of 
fitting functions and approximations to the full models dis- 
cussed here - however, we have tested extensively that the 
approximations and fitting functions used yield much smaller 



differences in the ultimate merger rates than the inherent un- 
certainties discussed here. 
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